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Abstract 

The creation and justification of the methods for minimax estimation of parameters of the exter- 
nal boundary value problems for the Helmholtz equation in unbounded domains are considered. 
When observations are distributed in subdomains, the determination of minimax estimates is 
reduced to the solution of integro-differential equations in bounded domains. When observations 
are distributed on a system of surfaces the problem is reduced to solving integral equations on 
an unclosed bounded surface which is a union of the boundary of the domain and this system 
of surfaces. Minimax estimation of the solutions to the boundary value problems from point 
observations is also studied. 

MSC-class: 35J25, 45Fxx, 45Kxx, 49xx 93E10 (Primary), 78M50 (Secondary) 

Introduction 

In the system analysis of complex processes described by partial differential equations 
(PDEs), an important problem is the optimal reconstruction (estimation) of parameters 
of the equations, like values of some functionals on their solutions or right-hand sides, 
from observations, which depend on the same solutions. 

These problems play an important role in various areas of science and engineer- 
ing. Depending on the character of a priori information, stochastic or deterministic 
approaches are possible. The choice is determined by the nature of the problem param- 
eters which can be random or not. Moreover, the optimality of estimations depends on 
a criterion with respect to which a given value is evaluated. 

The field of optimal control of PDEs has been strongly influenced by the work of 
J.L. Lions, who started the systematic study of optimal control problems for PDEs in 
[1], in particular, singular perturbation problems in [2] and ill-posed problems in [3]. A 
possible direction of research in this field consists in extending results from the finite- 
dimensional case such as Pontryagin's principle, second-order conditions, structure of 
bang-bang controls, singular arcs and so on. On the other hand partial differential 
equations have specific features such as finiteness of propagation for hyperbolic systems, 
or the smoothing effect of parabolic systems, so that they may present qualitative 
properties that are deeply different from the ones in the finite-dimensional case. The 
present study is devoted to a class of problems of optimal control and estimation for a 
specific family of PDEs of mathematical physics. 
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In practice, the data of boundary value problems (EVPs) for differential equations 
that simulate a physical or technological object are always given with uncertainty. For 
example, the right-hand sides of the equations, initial or boundary conditions may be 
known approximately; that is, they belong to certain bounded sets in the corresponding 
functional spaces. 

For solving the estimation problems we must have supplementary data (observations) 

y = C(f + r], 

where C is an operator that specifies the method of measuring and rj is the measurement 
error. As a rule, this error is not known and belongs to a certain given set and the 
operator is not invertible. Therefore, in general, from given it is not possible to 
uniquely reconstruct the sought-for solution (/? of a EVP and, consequently, quantity 
/((/?), where / is a given linear continuous functional. We see that a natural problem 
arises: to determine a quantity l{(p) which would provide the best (in a certain sense) 
approximation to the sought-for 

Let us briefly characterize the minimax approach to the solution of this problem. 
We are looking for linear with respect to observations optimal estimates of functionals 
of solutions and right-hand sides of EVPs based upon the condition of minimum of the 
maximal mean square error of estimation taken over the subsets mentioned above. 

These estimates were called minimax a priori or minimax program estimates (see 
[27], [28]). 

The situation when the unknown parameters of equations and observations are per- 
turbed by noise whose statistical characteristics are not known completely constitutes 
the case of special interest. 

In the absence of true information about distribution of random perturbations, the 
minimax approach proved to be a useful solution technique. This approach initiated and 
developed by N.N. Krasovskii [27], A.E. Kurzhanskii [[2$], O.G. Nakonechnyi [1], N.F. 
Kirichenko [TU], and E.M. Pshenichnyi enabled one to find optimal estimates of the EVP 
parameters for ordinary differential equations corresponding to the worst realizations 
of random perturbations. 

The present work is devoted to the creation and rigorous justification of constructive 
minimax estimation methods of parameters of the external EVPs for the Helmholtz 
equation in arbitrary unbounded domains with finite boundaries. We reduce the deter- 
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mination of minimax estimates to the solution of certain integro-diffcrcntial equations 
in bounded domains when observations are distributed in subdomains. When observa- 
tions are distributed on a system of surfaces (that simulate e.g. antennas) the problem 
is reduced to solving some integral equations on an unclosed bounded surface which is 
a union of the boundary of the domain and this system of surfaces. 

These estimation problem are of tremendous significance in many areas of applied 
electromagnetics, acoustics, contact mechanics. Therefore, comprehensive theoretical 
analysis of estimation techniques is an urgent task. 

Methods and objectives. The study is aimed at elaboration of the methods of 
guaranteed estimation of the values of linear functionals defined on solutions to external 
EVPs for the Helmholtz equation and their right-hand sides. 

This task can be fulfilled if the following problems are solved: 

• To reduce estimation of the values of functionals defined on the solutions to external 
EVPs and the right-hand sides of equations that enter the problem statement 
to certain problems of optimal control of systems that are described by certain 
conjugate EVPs for the Helmholtz equation in bounded domains with a quadratic 
quality criterion. 

• To obtain, for given restrictions on the unknown second moments of observation 
noise and unknown deterministic data of the EVPs under study, the systems of 
intcgro-differential and integral equations such that the minimax estimates of func- 
tionals are expressed in terms of their solutions 

• To prove unique solvability of the obtained systems of integro-differential and in- 
tegral equations. 

The object of study is observation problems under uncertainty when the functions 
that are observed on a system of subdomains or surfaces are coupled with the solutions 
to the considered EVPs via linear operators with additive measurement errors. 

The method of study. The systems of integro-differential and integral equations 
obtained in this work whose solutions are used to express minimax estimates arc based 
on the theory of generalized solutions to EVPs for the Helmholtz equation, utilization 
of the so-called Dirichlet-to-Neumann (DtN) data-transforming operators, and the the 
theory of potential in Sobolev spaces. 
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A remark on novelty. For the first time we consider tlie statement of tiie problem 
of minimax estimation of the parameters of external BVPs for the Helmholtz equa- 
tion with general boundary conditions that arise in the mathematical theory of wave 
diffraction. 

For the systems described by such BVPs, we obtain representations for minimax 
estimates of the values of functionals from the observed solutions and right-hand sides 
that enter the problem statement; quadratic restrictions are imposed on unknown deter- 
ministic data and second moment of observation noise. We also obtain representations 
for the estimation errors. The representations are obtained in terms of the solutions 
to certain uniquely solvable systems of integro-differential and integral equations in 
bounded domains. 

When the unknown solutions of the system states are observed that are described by 
external BVPs for the Helmholtz equation on a system of surfaces, we obtain systems of 
integro-differential equations in unbounded domains; the required minimax estimates 
are expressed via the solutions to these systems using integral operators of the potential 
theory in Sobolev spaces; and the BVPs are reduced to equivalent integral equation 
systems on multi-connected surfaces (or contours), the latter being a union of the 
obstacle boundary and the surfaces on which the observations are made. 

We prove the unique solvability of the obtained integral equations for any values of 
the wave number k such that Imk > 0, k ^ 0. 

Practical importance. The estimation techniques elaborated in this work are of 
big importance for the development of the theory of inverse acoustic and electromagnetic 
wave scattering by bounded obstacles. 

The methods and results of this study may be used for estimating under uncertain 
conditions of the system states described by BVPs for Helmholtz equation in more 
complicated domains with the boundaries that stretch to infinity (for example, in a 
domain K \ Q where is a layer between two parallel planes and is a bounded 
domain). In general, the developed estimation methods can be applied to obtaining 
minimax estimates of parameters for a wide class of problems of mathematical physics. 
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PART 1 



Minimax estimation of the solutions to the Helmholtz problems 
from observations distributed in subdomains 



1.1 Notations and definitions 

Let us introduce the notations and definitions that will be used in this work. 

) denotes a spatial variable that is varied in an open domain D C M"; 
is a Lebesgue measure in R"; 
x(M) is a characteristic function of the set M C W^] 
H^(W^) is a Sobolev space of index s: 

H'iW") = {ue L\W') : (1 + \y\Y^J^u{y) G L\W)}, 

where s > 0, ^^(M") is a space of square integrable functions in and Tu{y) denotes 
the Fourier transform of function u{x). If s < 0, then H~^{M?) denotes the space dual 
to H^iW'). Let D be a domain in (not necessarily bounded) with the Lipschitz 
boundary dD. Then ddD denotes the element of measure on contour dD. L^{dD) is a 
space of square integrable functions on dD] the function space 

LL{D) = {ue V\D) : u\Dmn e L\D n ^Ir) 

for every > such that D n ri/j ^ 0} , 

Introduce also the Sobolev spaces with the corresponding norms: 

H\D) = {u\d : u e H'{R^)} (s G R), 



u 



H\dD) = < 



dD-. ue H'+'/^R^) (s > 0), 

L\dD) (s = 0), , 

H~-^{dDy (dual space with respect to H~-'{r)) (s < 0) 



for every R> such that D f) Qr ^ , 
H\D, A) := {u G D'{D), u G H^[D), Au G L\D)}, (LI) 
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HUD. A) := {u e D\D), u G HI{D). Au e LI{D)}, (1.2) 

u is identically zero outside some ball centered at the origin} , (1.3) 

where V'{D) is the space of distributions in D; here and below by ^Ir we denote the 
ball flu := {x : \x\ < R}; the Laplacian is taken in the sense of distributions in D; and 
s G R. 

Theorem (The trace theorem for H^{D), [5], p. 102). 

For any Lipschitz domain D an operator 7/5 : C^{D) C^{r) can be extended to 
a continuous and surjective operator '■ ^iL(^) ~^ H^^^{dD). 
We denote by 77V the Neumann trace operator 

{'^nu){x) (grad w(a:), n(a:))Mn, x G dD, u G C^{D), 

: C^{D) C^\dD) and can be extended to a continuous and surjective operator 
77V : Hl^^{D, A) H^^I'^{dD). This operator will further be denoted by d/dv. 

< •, • > H-^i'^{dD)y.m/'^{dD) denotes the duality relation between spaces H^^^'^{dD) 
and H^/\dD), which is an extension of the inner product in L^{dD) in the following 
sense: if r G L'^{dD), then the following relation holds 

<r,w >H~i/^{dD)xm/^{dD)= / rwddD \/w ^ H^/'^{dD). * 

JdD 

Let be a Hilbert space over the set of complex numbers C with the inner product 
(•, ■)h and norm || • \\h. By L'^iTi^ H) we denote the Bochner space composed of randomS 
elements ^ = ^(cj) defined on a certain probability space (E,;B, P) with values in H 
such that 

ii^iiW)= / umidP{u)<^. (1.4) 

In this case there exists the Bochner integral := dP{u) G H which is called 

the mathematical expectation or the mean value of random element ^{uj) and satisfies 
the condition 

{hM)H = J {h,^{u))HdP{u) yheH. (1.5) 

^Random element ^ with values in Hilbert space H is considered as a function ^ : E ^ iJ imaging random events 
E B to Borel sets in H (Borel cr-algebra in H is generated by open sets in H). 
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Being applied to random variable ^ this expression leads to a usual definition (value) of 
its mathematical expectation because the Bochner integral ( 11.41 ) reduces to a Lebesgue 
integral with probability measure dP{uj). 

In L^(E, H) one can introduce the inner product 

Applying the sign of mathematical expectation, one can write relationships (11.41) — ( fL6l ) 
as 

ll^llW) = E||^(a;)||L (1.7) 

(/i,E0i7 = E(/i,^(a;))i/ V/i G (1.8) 

(^, T/)L^(E,i7) := E(^(a;), t^{uj))h V^, t/ G L2(E, H). (1.9) 

LP'{Ti^H) equipped with norm (11.71) and inner product (11.91) is a Hilbert space. 

Consider the problem of finding a solution to the exterior Neumann problem for the 
Helmholtz equation. 

Assume that Q G is a bounded domain such that = F is a Lipschitz contour 
and r^o is a bounded subdomain of \ Q, C \ Cl. Given a function / defined in 
the domain \ H such as / = outside Hq, / ^ ^^(^o) and a function g G i7^^/^(r), 
find iIj G HlX{R^ \ ^), A), such that 

- (A + A;2)V;(x) = /(x) in (1.10) 

dil) 



du 



= g on (1.11 



= o(l/r^/2), r = \xl r ^ oo (1.12) 



dr 

with an equivalent variation formulation: find ifj G H^^^i^ \ Cl) such that 



{V^Ve -k'^e)dx= / fedx+ gOdV (1.13) 

for all 9 G H^^^^(M?\Q) such that i/j satisfies the Sommerfeld radiation condition (11.121) . 
Here we suppose that k is the wave number with Im /c > 0, /c 7^ 0. 

Consider also the following problem: find i/j G H^{{Qji \ Q), A) such that 

- (A + e)ilj{x) = f{x) in Qr \ n, (1.14) 
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^ = 5'onr, (1.15) 

^ = M^^V on r^, (1.16) 

where / G L2(Qo) and g G //-^/^(r), : H^''^{Tr) H-^/^Tr) is the Dirichlet- 
to-Neumann map (DtN map) defined by 

{Mi'WR,e) := A^^|!M r^(i?,0)e-(^-^)rf0, (1.17) 

and is a large disk containing Q and supp /. It is known that problems (11.101 ) — ( [1.121) 
and (11.141) — ( Tl. 161) are equivalent in the following sense (see [6] — [H]). If ^/^ G Hl^^{{M.'^ \ 
Q), A), is a solution of (11.101) — ( [1.121) . then the restriction of i/j to ^Ir \ Q belongs to 
H\{nR\n),A) and is a solution to (fLTID - flTlB . Conversely, ifi/jE //i((Qi?\ H), A) 
is a solution to (11.141) — ( [1.161) , then this solution extended to the domain M? \ CIr by 

i^{rp,ep) = J^IJM rv^o(i?,0)e-(^-^)rf0, TP > R, (1.18) 

belongs to Hl^{{R^\n), A) and satisfies flLTm - flTl^ . Here 0) ■= ^{R^ 0) is the 

trace of the solution to problem (11.141) — ( [1.161) on Tr and (rp, 6p) are polar coordinates 
of the point P G \ Qr. 

To formulate an equivalent variational setting of problem (11.141) — ( [1.161) , we introduce 
the continuous sesquilinear form a(-, •) : H^{flR \ Q) x H^{flR \ H) ^ C defined as 

r(l) 

ivi 

R 

Then an equivalent variational formulation of problem (11.141) — ( jl. 161) can be written as 
follows: find ip G H^{Qr \ Q) such that 

a{iij,e) = i{e) yeeH\nR\n), (1.20) 

where 

1(9) := / fedx+ f gOdV (1.21) 

is a continuous semihnear functional on H^{Q.r \ Q). 

In order to set an adjoint problem of (11.141) — ( [1.161 ) which will be used below we 
introduce a sesquilinear form 

a*(v^,^):= / {Vijve-Pije)dx- [ Mf^iijedVp, (1.22) 



a{i),e)-- {ViljV9-k^iP0)dx- Mlr't/jedTR. (1.19) 
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k 



where Aff * : H^'^iTR) H-^/^{rie) is the map defined by 

(MfM(i?,(^) := l^^gl r«(ij,0)e»<*-«d^. (1.23) 

Lemma 1.1. The sesquilinear form a*{il^,6) is an adjoint ofa(il^,9). 
Proof. Defining 

JflR\Cl 

and 

we liave 

a(V^, 6) = ai(V^, 6) - a2(V^, 9). (1.24) 

Obviously, 



a*(^,^) :=ai(^,V^) = / [VOV^j - m ^j) dx = / {Vt/jVO - k'^jj 6) dx. (1.25) 
Taking into account that 



we find 

4(V^,^) :=^^(M0 = / Ml'^ei^dTR= f Ml'^eil^dTn 



Jo m^HkR) Jo 



= R r^(i?,a)e-(>^-)rfaV^(i?,x)^X 



'277 jT o-(2)VjUD\ /•27r 



= ^ r^T. \S^! t\iR,x)e'''^''-''^dxe(R^da= [ Mf^^OdTn. 

n€Z Hn {kR) Jo JTr 

(1.26) 
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From ( Oil )-( OBl ) and the equality 



□ 



we obtain the required assertion. 

Now we can state the variational problem adjoint of (11.201) : 

Given functions f and g introduced on pagelB, find i/j G \ Q) such that 

a*{il),6) = l{6) W e H\nR\n), (1.27) 

where functional l{9) is defined by (11.211) . 

It is easy to see that variational problem (11.271) as well as the following problems: 
(i) find V^G//i((^i?\^),^) such that 

- (A + P)V^(x) = f{x) in Qr \ (1.28) 

g on r, (1.29) 



dv ^ 

and (ii) find G //iL((M^ \ f^), A) such that 



- Mf V on r^, (1.30) 



(A + k^)ilj{x) = f{x) in \ Q, (1.31) 
=g on r, (1.32) 



+ iki/j = o{l/r^^'^), r = |a;|, r oo, ecjin Im/c > (1.33) 



dr 

are equivalent and for any R there exists a positive constant a > independent of / 
and g (but dependent on R) such that 

< oi{\\f\\H-i{n,) + Iklli7-i/2(r))- (1-34) 

If any of the data in EVPs flLll) - f|LT6l) or ( OSl ) - ( OOl ) is random (e.g. forcing 
function / or Neumann boundary data g), then the solution i/j will be a random function 
and corresponding stochastic EVPs are formulated as follows. 

Given / G L2(E, L^Qq)), g G L2(E, //-i/2(r)), find ^ G L2(E, H^Qr)) such that 

Ea(^, 6) = El{e) W G L2(E, H\nR \ H)), (1.35) 
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and 
where 

E/(^):=e| / f{x,u)e{x,Lu)dx^ f g{-,u)e{-,Lu) dvl (1.37) 

is a continuous semilinear functional on L'^{H^{flji \ Q)), Ea(-, •) and Ea*(-, •) : 

L2(E, L'^{H\nR \ Q)) X L2(E, L'^{H\nR \ Q)) ^ C are continuous sesquilinear forms 
defined as 



Ea(?/;, ^) := E <^ / (V^(a;, a;) V^(a;, cj) - k'ij{x, u) 6{x, u) dx 



^ M(^V(-,cu)^(-,a;)dr^| (1.38) 



and 



Ea*(V^,^) := E<^ / {Vilj{x,u)Ve{x,u) - k^ilj{x,u) e{x,u) dx 



^ MfV(-,a;)^(-,a;)rfr^|. (1.39) 



It is known that problems (11.351 ) and (11.361) have unique solutions and there exists 
a positive constant a > independent of / and g such that 

II^IU2(s,i7i(r!fl)) < «(||/||l2(s,l2(Oo)) + ||5'llL2(i],i/-i/2(r)))- (1-40) 



Such problems were investigated in [12] and [TH], including the construction of finite 
element methods of their numerical solution. 

Problem (11.351) is equivalent to the following ones: find ijj G i^^(S, H^{{^}r \ Q), A)) 
such that 

- {A + k'^)ij{x,uj) = f{x , to) in Qi? \ (1-41) 



du 



= g{-,iu) onT, (1.42) 



?M = MlV,.)onr„ (1.43) 
or find G L2(S, H^Qr \ Q)) satisfying 

a{^p{-,uj),e) = i{e) ye g H\nR\n). (1.44) 
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Analogously, problem ( 11.361) is equivalent to the following problems: find ijj G H^{{^r\ 
H), A)) such that 

- (A + P)7/'(x,a;) = f{x 
5^(-,w) 



5z/ 



^(•,a;)onr, (1.46) 



M^ = MfV(,c.)onr„ (1.47) 
or find ^ G L2(E, //^(Qi? \ Q)) satisfying 

a*(^(-,a;),^) = /(^) G //^(Qi? \ Q). (1.48) 

The right-hand sides in (11.411) — ( [1.481) are considered for every realization of random 
fields f{-,uj) and g{-,u) which belong with probability 1 to the spaces L'^{Qq) and 
L^(r), respectively, and the equalities are satisfied almost certainly. 

1.2 Statement of the estimation problem 

Consider the exterior Neumann problem for the Helmholtz equation: find a distribution 
if G I>'(R2 \ Q) such that 

^G//L((M^Q),A), (1.49) 
- {A + k'^)ip{x) = f{x) in R^\n, (1.50) 
^ = ^onr, (1.51) 

dip 

— ikif = o(l/r"^/^), r = \x\, r ^ oo, Imk > 0, (1.52) 

or 

where k is the wave number with Im /c > 0, / is a source term distributed in bounded 
subdomairti in \ / ^ L^{no), and g G L^{T). 

BVP (11.491) — ( [1.521) simulates, in particular, acoustic or electromagnetic scattering 
from an infinite sound-hard (perfectly conducting) cylinder with cross-section fl. 

Denote by Go the set of pairs of functions (/,^) satisfying the inequality 



/ Qi{f-fo){x){f-fo){x)dx^ Q2{~9-go){~g-9o)dr<i, (1.53) 

and by Gi the set of random functions ^(•) = (^i(-), • • • ,^m(")) defined on ili x • • • x 
with integrable second moments E|^/fc(x)|^ satisfying conditions 



E^k{x) = 0, k = l,m. (1.54) 



^This means that the function / is defined in the domain ]S? \ Q, f = outside flo, and / G L^(rio) 
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lit p 

J2 / B\ik{^)\'rl{x)dx<l, (1.55) 
k=i "^^fe 

where /o is defined in the domain R^\n, /o = outside Hq, /o £ ^^(^o) and 5^0 ^ ^^(r), 
/o and are prescribed functions, Qi and Q2 are Hermitian operators in L^(Qo) and 
L^(r), respectively, for which there exist bounded inverse operators Q^^ and Q2 
rk{x): /c = l,m, are nonvanishing functions continuous on sets Qk- 

We suppose that functions f{x) and ^'(x) in equations ( 11.501 ) and ( 11.511 ) are not 
known exactly; it is known only that (f^g) G Gq. 



Assume that in subdomains Qk, k = l,m, of domain Q the following functions are 
observed 



Vk 



{x) = gk{x,y)(p{y) dy + ^k{x,u), x e Qk, k = l,m, (1.56) 

where (p is a solution of BVP ( 11.491 ) — ( [1.521 ), gk{x,y) G L'^{^k x ^fc) are prescribed 
functions and ^fc(a:,ct;) are the choice functions of random fields ^fc(x) with unknown 
second moments such that ^{■) = (^i(-)? • • • ? ^m(")) ^ ^i- 

Let Iq be a given function defined in a bounded subdomain uq C Q belonging to 

The estimation problem consists in the following. From observations ( 11.561 ) of the 
state <-p{x) of the system described by BVP ( 11.491 ) — ( [1.521 ) under conditions ( 11.541 ) — ( [1.551 ) 
it is necessary to estimate the value of the linear functional 

l{(p) = / loix)(f{x) dx (1.57) 

JUJQ 

in the class of the estimates linear with respect to observations which have the form 

K^) = ^ Uk{x)yk{x) dx + c, (1.58) 
k=i -^^^ 

where Uk G L^(Qfc), k = l,m, c G C. 

Denote by w = (wi, . . . , Um) an element belonging to H := L'^{Qi) x . . . x L^(Qfc). 

Definition 1.1. An estimate l{(p) is called a minimax estimate of the l{(p) if an element 
u = (wi, . . . , Um) G H and a number c G C are determined from the condition 

sup nm-mf^ mf 
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Here (p is a solution to problem (l{1.4S\ )— (il.5^) when f = fj g = and l{(p) 



EfcLi Uk{x)yk{x) dx + c, yk{x) = J^^ gk{x, y)(p{y) dy + ^k{x), x e Qk: k = 1, m, 
The quantity 

a:={ sup E[/((^) - (;^_59) 

is called the error of the minimax estimation ofl{(p). 

Thus, the minimax estimate is an estimate minimizing the maximal mean-square 
estimation error calculated for the "worst" implementation of perturbations. 

1.3 Reduction of the estimation problem to an optimal control problem 

Lemma 1.2. The problem of finding the minimax estimate ofl{(p) is equivalent to the 
problem of optimal control of a system described by a EVP 

z{-,u)eH\{nR\n),A), (1.60) 

m ^ 

A + k^)z{x;u) = Xujo{^)k{x)-y2xnk{x) / gk{r],x)uk{r]) dr] in UrXQ, (1.61) 

OonT, (1.62) 



k=i 
dz{-] u) 



dv 

dr 

with the quality criterion 



hrfz(--u)onTR (1.63) 



/(«):= /QrM.r«)»<ix+/Q-.(.;«);0^<ir 

m „ 

+ y / rf{x)\uk{x)\^dx inf, (1.64) 
k=i ^^^fe 



where R is chosen so that Q^i C \ i = 0,m, uq G \ 
Proof. Taking into account (11.561 ). (11.571) . and (11.581) . we obtain 

p m „ 

l{(p) - l{(p) = I lo{x)(p{x) dx - ^ Uk{x)yk{x) dx - c 
J ujQ J 

= / l^{x)(p{x) dx -^2Z / ^^^^^ / 9k{x,y)(p{y) dydx / Uk{x)^k{x) dx - c 
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= / lo{x)(p{x) dx / / gk{y,x)uk{y)dyip{x)dx / Uk{x)ik{x) dx - c 

Jloq j^^^ J ilk Jrik i^^Y '^^k 

' m p \ 

Xuo{x)k{x) / 9k{y,x)uk{y)dy\^{x)dx 
'nit\n \ J ilk I 

m „ 

~ XI / '^k{x)ik{x) dx - c. (1.65) 

For any fixed u = (wi, . . . , G iif introduce the function z{x] u) as a unique solution 
of problem ( 11.601 ) — ( [1.631 ) . According to the equivalent variational formulation of this 
problem, it means that z{x] u) satisfies the integral identity 

e)= / {Vz{x]u)Ve{x)-k^z{x;u)e{x))dx- Mfh{-]u)edTR 
JnR\n Jtr 

= \ Xioo{x)lo{x) -y2^^ki^) gk{y,x)uk{y)dy\ e{x)dx ^6 e H\nji\n). 

(1.66) 

Set 6 = (p m ( II ■661 ). Then we obtain 
a-(.(.;.),^)= / (V.(.x;«)V^-P.(x;.)vKi)dx- / Mf .(.;«)^<ir« 



m 



Xujo{x)lo{x) -y^^^ki^) 9k{y,x)uk{y)dy] (fdx. (1.67) 

On the other hand, since is a solution of problem ( 11.491 ) — ( [1.521 ) with / = / and 
g = setting i/j = (p and 6 = z{-; u) in ( 11.191 ), we find 

JnR\n Jtr 



= / f{x)z{x]u) dx + I gz{-;u)dT. (1.68 
Jrio Jr 



By Lemma 1.1, a*{z^(p) = a{(p^z). This identity and ( 11.641 ), ( 11.671 ), and ( 11.681 ) imply 



l{(f)-l{(p)= \Xujo{^)k{x) -y2^^k{^') 9k{y,x)uk{y)dy\(p{x)dx 

= / f {x) z{x; u) dx + / gz{-]u) dV — '^^ / Uk{x)^k{x) dx — c. (1.69) 
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Taking into consideration the relationship Dr/ = E|?7 — Er/p = E|?7p — |E?7p that 
couples dispersion Dij of the complex random variable t] = t]i + ir]2 and its expectation 
Er/ = E?7i + iE?72, we obtain from the last formulas 

m 2 



E 



/ f [x] z{x; u) dx + I gz{-]u)dr 
'rio Jr 



+ E 



/ Uk{x)ik{x)dx 
k=i '^^'^ 



Therefore, 



inf sup E\l{(f) - = 



= inf sup 

"^'^ {f,g)eGo 



/ f {x) z{x; u) dx + I gz[-;u)dT — c 
'nn Jr 



sup E 



m « 

/ Uk{x)ik{x)dx 



(1.70) 



In order to calculate the first term on the right-hand side of ( 11.701) make use of the 
generalized Cauchy— Bunyakovsky inequality in ( 11.531) . We have 

/. „ 2 



inf sup 



/ f {x) z{x] u) dx + I gz{-;u)dT — c 



= inf sup 

(/,.g)eGo 



< 



{f{x) - fo{x))z{x]u) dx + / {g - go)z{-;u)dr 

+ / fo{x)z{x)dx+ I gQz{-]u)dT — c 



X <^ / Qlif - fo){x){f{x)- fo{x))dx^ / Q2{g-gQ){g-go)dr\ 
IJno Jr ) 

< / Qi^z{x;u)z{x;u) dx + / Q2^z{-;u)z{-;u) dT. (1.71) 

The direct substitution shows that that inequality ( 11.711) is transformed to an equality 
on the element {f''^\-):g^^^): where 

/W(x) :=^Q^h{x;u) + Mx), 
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~9^'^ ■.= \Q^'z{--u)+g,, 



d 



1/2 



Therefore 



inf sup 



{f{x) - fo{x))z{x; u)dx+ {g - go)z{-; u) dT - c 



= / Qi^z{x;u)z{x;u) dx + / Q2^z{-;u)z{-;u) dT (1-72) 
Jn Jt 



with 



c= z{x] u)fo{x) dx / z{-;u)godr. 

In order to calculate the second term on the right-hand side of ( 11.701 ). note that the 
Cauchy— Bunyakovsky inequality and ( 11.551) yield 



X] / Uk{x)ik{x)dx 



< 



m 



/ rf^^{x)\uk{x)fdxY^ / rl{x)\ik{x)\'^dx, 

, T J ill. 1, 1 J ilk 



k=l 



k=l 



the latter implies 



sup E 

led 



m p 

/ Uk{x)lk{x)dx 



m 



r^2(a:)|wA:(x)p 



fc=i ^^'^ 



However 



E 



o 

m „ ^ m „ 

J2 Uk{x)if\x)dx =J2 ^kHx)\Mx)f 

k=l ^^'^ yfc=l ^^fe 



where 



iyr^'^{x)uk{x] 



x G ilfc, k = l,m, 



|<»' = {ll\),...,|(?)(.))eGi, 
and z/ is a random variable with Ez/ = and E|z/|^ = 1. Therefore 

^ 2 ^ 



sup E 

led 



110 n 

/ Uk{x)ik{x)dx 
k=i -^^^ 



= E / (1-73) 



which proves the required assertion. The validity of Lemma 2.2 follows now from 
relationships ( 11701 ). ( fTT^ . and ( IITHl ). □ 
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1.4 Representation of minimax estimates and estimation errors 

In the course of the proof of Theorem 1.1 below, we show that the solution to the 
optimal control problem ( 11.21 ) — ( IT. 641 ) (and therefore, the determination of the minimax 
estimate, in line with Theorem 1.2) is reduced to the solution of a certain integro- 
differential equation system. Namely, the following statement holds. 

Theorem 1.1. The minimax estimate of l{(p) has the form 

. . m p 

^('^) = X] / Uk{x)yk{x) dx + c, (1.74) 

where 



Uk{x) = rl{x 



/ 9k{x,y)p{y) dy, k = l,m, c= / z{x)fo{x) dx+ / zgodT, (1.75) 

and functions z and p are determined from the solution to the following problem: 

zeH\{nR\n),A), (1.76) 

m „ 

- {A + k^)z{x) = Xujo{x)lo{x) -^X^ki^) gk{r],x)uk{r])dr] innR\n, (1.77) 

J fij. 

OonT, (1.78) 



k=i 
dz 
du 



^ = Mf^z on r^, (1.79) 

l?G//i((Qi?\Q),A), (1.80) 

- (A + k')p{x) = xn,{x)Q^h{x) m \ Q, (1.81) 

^ = Q^^z, onT, (1.82) 

^ = mI'^p on Fr. (1.83) 

Problem (11.761 ) — ( [1.831 ) is uniquely solvable. The restrictions of the solutions of this 
problem corresponding to R = Ri and R = R2 on ^i?=min{i?i,i?2} coincide with the 
solution corresponding to R = min{i?i, R2}. 

The estimation error a is determined by the formul(^ 

a = l{py/^=(^J Iq{x)p{x) dx^ . (1.84) 



■^In the proof of this theorem, we show that the value of J^^ lo(x)p(x) dx is real. 
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Proof. Let us show first that the optimal control problem ( 11.601) — ( [1.641) is uniquely 
solvable; that is, there exists one and only one element w G i7, at which functional 
( 11.641) attains the minimum value, I{u) = infuen I{u). 

It is easy to see that the solution z{x; u) of BVP ( 11.601) — ( [1.631) can be represented 

as 

z{x]u) = zq{x) + z{x]u)^ (1.85) 
where zq{x) and z{x] u) are solutions of the following EVPs 



- (A + k^)zQ{x) = Xu{x)Iq{x) in VLr \ 

— — = on r, 

ov 

-— = Ml 'zq on Tr 

and 

z{--u)eH\{^R\Q),/\), 

m « 

- {A + k^)z{x]u) = - ^Xnui^) / gk{r],x)uk{r]) df] in nji\n, 



dz{-]u) 

— - — = on r, 



1.86) 

1.87) 
1.88) 

1.89) 

1.90) 
1.91) 

1.92) 
1.93) 



— — — = Ml 'z[-] u) on Tr. 
Using representations ( 11.851) for z{t] w), write functional I{u) in the form 

I{u) = I Qi^ z{x]u)z{x]u) dx + I z{-]u)z{-]u) dV + '^^ I r^'^{x)\uk{x)\^ dx 
J Qq Jr J rik 

= / Qi^{zo{x) + z{x;u)){zo{x) + z{x;u))dx + / Q2^{zo + z{-;u)){zo + z{-;u)) dT 

m „ 

/ rj^'^ix)\uk{x)\'^ dx = i{u) + L{u) + Co, 



where 



I{u) := / Qi^z{x] u)z{x] u) dx / Q2^^5(-; 'w)i(-; w) c/F + / 



rf^'^{x)\uk{x)\'^ dx 
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is a quadratic functional in the space H which corresponds to a semi-bihnear continuous 
Hermitian form 

Jrio 

m „ 

+ / r^'^{x)uk{x)vk{x) dx (1.94) 
fc=i ^^'^ 

on H X H and satisfies 

I{u) > a\\u\\'jj^ \fu E H, a = const; (1.95) 

note that 

L(w) := 2Re / Qi^z{x] u)zo{x) dx 2Re / Q2^z{-]u)'zodr 
Jrio Jr 

is a Hnear continuous functional in H and 

Co := / Qi^zo{x)zq{x) dx ^ / Q2^zo'z^dr. 

Prove, for example, the continuity of form (11.941) : namely, the inequality 

|7r(f , w)! < cll-ulli/JIwllij^ V-UjWGV, c = const (1.96) 

(the continuity of linear functional L{v) is proved in a similar manner). 

Using estimate (11.341) obtained above and the Cauchy— Bunyakovsky inequality we 
have 

\7r{u,v)\< I / Qi^z{x;u)z{x;u) dx \ I / Qi^z{x;v)z{x;v) dx \ 
\Jno J \Jno J 

+ (^J^Q2'z{--u)J{^dr^ ' {jQ2^z{--v)l{^)dv\j ' 

< / Qi^z{x;u)z{x;u) dx + / Q2^z{-;u)z{-;u) dT + 2. / r'i^'^{x)\uk{x)\'^ dx \ 
\Jn„ Jt Jnk J 

^1/2 

X I / Qi^z{x;v)z{x;v) dx + I Q2^z{-;v)z{-;v) dT + y I r^'^{x)\vk{x)f' dx 



ft n I I L n 

/ Qi^z{x]v)z{x]v)dx+ / Q2^z{-]v)z{-]v) dV + '^^ / 
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\ 1/2 



\Q^^z{x;u)fdx] ' ( / \z{x]u)fdx] +( f \Q2^z{-,u)f dv] 
'no J \Jno J \Jt J 

1/2 / ^ „ \ l/2>, 1/2 



\ 1/2 



< 



X 1/2 / rn \ / m \ 



X 1/2 



|Qr^^(^;^)Pf^2:) ' ( / \z{x]v)fdx] +( f \Q2^z{-]v)fdr 

\l/2fm^ \ 1/2 / m „ \ 1/2 ^ 1/2 



1/2 



X 



/ \z{x;u)\^dx+ I |5(-;w)p(ir + V / 



1/2 



1/2 



x< / \z(x;v, 
'no 



n III n 

/ \z{--v)fdV^^ \ \vk{x)fdx 

1 /9 

< max{||Q^^||, ||Q^^||,/?}{P(-; + Uto^Is w)|i2(r) + ||w|||^} 

X {ll^(-;^)||?/i(f7o) + ll7o5(-;^)li2(r) + 
where ^ := maxl<fc<rr^ max^^f^^ ^fc^(a^) > 0. 



(1.97) 



Setting in (|1.34|) = / = -YA^^iXn^i^) J^^gk{r],x)ukir]) drj, and 5^ = 0, 

we find 



< a 



^Xn^ix) / gkir],x)ukir])dr] 
1 1 <^ f^fc 



A:=l 



<ci||w||,f, (1.98) 



where ci is a constant independent of i?. The trace theorem and inequahty (11.981 ) imply 

\\1oz{-]u)\\l2(t) < C2\\z{-;u)\\Hi{nn\n) < C3\\u\\h, 02,03 = const. (1.99) 
From I^M)-I^M) we have 

\7t{u,v)\ < max{||Q^^||, 11^2^11' /^KciPIll/ + C4||w|||f + ||w|||-)^^^ 

X (cilbllff + OiWvWjj + plllf)^/^ < c||w||i/p||i/, 
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where c is a constant independent of R. 

Thus inequahty ( 11.961 ) and, consequently, the continuity of form ( 11.941 ) are proved. 

In Hne with Remark 1.4 to Theorem 1.1 proved in |1], p. 11, the latter statements 
imply the existence of the unique element u ^ H such that 

inf /(^). 

Therefore, for any r E R and v E H, the following relations are valid 



-^I(u + Tv) 
ar 



= and —liu + iTV. 
r=o dr 



= 0, 



r=0 



(1.100) 



where i = Since z{x] u + tv) = z{x; u) + tz{x] v), where z{x; v) is the unique 

solution to BVP ( ILM ) - ( fLM ) at u = v and Iq = 0, the first relation in ( frTOOj ) yields 



= lim — I 



(Qi ^z{-; u + tv), z{-] u + Tv))L^{na) - {Q^ ^z{-; u), z{-; u))L^{no) 
+ u + tv), z{-; u + tv))l2[t)-{Q2^^{-'^ ^), z{-; u))l^t) 

m 

m 



Similarly, taking into account that z{x; u + itv) = w) + iTz{x; v), we find 

1 (i ,^ M 

= -—I{u + lTV)\r=0 

2 CLT 



= Im|(Q/z(-; w),5(-;?;))l2(Oo) + {Q2^z{-;u), z{-;v))t + ^ (^fc ^^fc' ^fc)L2(f],)}; 

i=l 

consequently, 

TO 

(Qr^^(-; u), z{-; v))L^n,) + (Q2 ^^(S u), z{-; v))t + ^ {rt^Uk, Vk) ^2^^^) = 0. (1.101) 
Introduce a function p{x) G H^{{flji \ Q), A) as the unique solution to the BVP 



(A + k^)p{x) = XOo(a^)Qi ^z{x] u) in VLr \ 



(1.102) 
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^ = Q2^z{-;u) onr, (1.103) 
^ = Mj^^p on r^, (1.104) 



or to an equivalent variational problem 



a{p,e)= I {Vpve-k'pe)dx- I M^^^pedVR 

JnR\n Jtr 

= 1 xn,{x)Qi^z{x-u)edx+ f Q^'z{--u)edTR yeeH\nR\n). (1.105) 

JnR\ri Jtr 



Setting in (11.1051) 9 = z{-;v), we obtain 



JnR\ri Jtr 

= / xno{x)Qi^z{x]u)z{x]v)dx ^ / Q2^z{-]u)z{-]v) dTji. (1.106) 
JnR\n Jtr 

Taking into account the fact that z{-;v) satisfies variation equation (11.271) with i/j = 
z{-; v) equivalent to BVP (ll.90l) -( fT.93l) with u = v and putting in (11.271) 6 = p, we have 

«•(!(.;.).,)= / (VI(,x;.)VM..)-P.(..;.)p(..))<i- / M^,(,.;v)pdr, 

JflR\Cl JTr 

p m p 

= -/ _^Xnk{x) gk{y,x)vk{y)dyp{x)dx. (1.107) 



Since a*(5(-; v),p) = a{p, z{-] v)), we have 

/ Qi^z{x]u)z{x;v) dx / Q2^z{-]u)z{-;v)dTji 



/ y^Xr!fe(a^) / gk{y,x)vk{y)dyp{x)dx 



= -^/ (/ gk{x,y)p{y)dy\vk{x)dx. (1.108) 



From ( frnm it follows that 



m p 

{Qi^z{-] u), z{-] v))L2{n^) + {Q2^z{-] u), z{-] v))t = - XI / 

k=l ^^fe 



r^'^{x)uk{x)vk{x) dx. 

(1.109) 
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Relations ( OHHD and ( OIM ) imply 

m p 

/ r^^{x)uk{x)vk{x) dx 
k=i -^^^ 

m r / r 

= E 

Hence, 



rf?/ ^^fc(3:) dx>Jvk G L (Q^:), k=l,m. 



k=i 



ilk \J / 



Uk[x) = rl{x) I gk{x,y)p[y)dy. 



Now let us establish the validity of formula (11.841) . We have 



:= sup M[l{(p) - l{(p)f = 
(/,5)eGo,|eGi 



/ Qi^z{x)z{x) dx ^ / Q2"^2Z(ir + ^^ / 



gk{x,y)p{y) dy 



Transform the sum of first two terms. Make use of equalities ( 11.1321) — ( 11.1351) to obtain 



a{p,z)= / xno{x)Qi^z{x)z{x) dx + / Q2^zzdr, (1.110) 



hence 



m „ 



fl{x] 



9k{x,y)p{y) dy 



dx. 



(1.111) 



Note that z satisfies ( 11.761) — ( Tl. 791) which yields an integral identity 



«o(x. 



m 

-I]xr^.(2:) / gk{y,x)uk{y) dy]e{^dx ^6 e H\nR\n); (1.112) 
setting 6 = p, we find 



,1 



a%z,p)= \ Xu;o{x)k{x) -y^^^ki^) gk{y,x)uk{y)dy\p{x)dx. 



From the latter relations, the formula a*{^z^p) = a(]7, z), and ( 11.1111) it follows that 
CT^ = / lo{x)p{x) dx - ^ / gk{y,x)uk{y) dyp{x) dx 
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k=i ^'^ 



dx. 



However, 



therefore, 



Uk[y) = rliy) / Qkiy, r])p{r]) dr], 



m 



5Z / / 9kiy,x)uk{y)dyp{x)dx 



k=l 



J ^k 



^ / 9kiy,x)rl{y) gk{y,r])p{j])dif]dyp{x)dx 

k=l ''^k ''^k J^k 



= X] / / 9kix,y)p{y)dy gkix,Tj)p{T]) drj dx 

k=l '= 



Finally, we obtain 



= / rl{x) / gk{x,y)p{y)dy 

k=l 

(7^ = / /o(2^)]5(x) C^X. 



□ 



In the following theorem we obtain an alternative representation for minimax esti- 
mate l{(p) that does not depend on the form of functional 



Theorem 1.2. The minimax estimate of l(ip) has the form 
where function is a solution to the following problem: 

peL\j:,H\inR\n),A)), 



(1.113) 



(1.114) 



{A + k^)p{x,uj) 

= Xnkix)Y^ / rl{T)gk{T,x) 
k=i '^^^ 



ykir,uj)- (f{ri,uj)gkiT,ri)dri 

'J ^k 



dp{-,u) 



d 



= onT, 



dr in Qr \ il, 
(1.115) 

(1.116) 
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(^GL2(S,//i((^i?\^),A)), 
(A + k^)^{x, uj) = Xno{x)Qi^p{x, u) + fo{x) in Qr \ 

= Q2 P{-,^)+9o, on r, 



Ml''(p{-,u) on Tr, 



^ ^(1) 



1.117) 
1.118) 
1.119) 

1.120) 
1.121) 



where equalities ({1.1151) — jil.llTl) and Iil.ll9\) — ^1.121\) are fulfilled with probability 1. 
Problem ( 11.1141) — ( 11.1211) is uniquely solvable. The restrictions of the solutions of this 



problem on ^^i?=min{i?i,i?2} corresponding to R = Ri and R = R2 coincide with the 
solution corresponding to R = min{i?i, i?2}- 

Proof. The proof is similar to the that of Theorem 2.1. Consider the problem of 
optimal control of the equation system 

p{;-,u)e L2(S, H\{nR \ Q), A)), (1.122) 
{A + k'^)p{x,u;u) = do{x,u;) / xn^i^) gk{r, x)u{T,u)dT in Qr\Q, (1.123) 

= on r, (1.124) 



dp{-,uj; u] 



dv 

dp{-,u;;u) (2) 

— = >(•, cj; u) on Tr (1.125) 



with the cost function 



I{u)=E< / Q^^{p{-,u;u) + Qifo){x){p{-,u]u) + Qifo){x)dx 
Uno 



+ J Q2^{p{-^^'^u) + Q2go){x){p{-,uj;u) + Q2go){x) dV 

+ / r.\x)\uk{x;cj)fdx\ inf . (1.126) 



where 



m „ 

do{x,u)) ■.= ^ / xnk{x)gk{r,x)rl{T)yk{T](p,^k{^))dT, 
k=i ^^'^ 
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The form of functional I{u) and proof of Theorem 2.1 suggest that there is one and 
only one element u G H) such that 

I(u) = inf I(u). 
Next, denoting by (p{t; to) the unique solution to the BVP 

-{A + k'^)(f{x,u) = xno{x)Qi^p{x,u]u) + fo{x) in Qr\Ci, 

Q2^pi-,uj;u) + go, on T, 



and making use of virtually the same reasoning that led to the proof of Theorem 2.1 
(by applying estimate (11.401) instead of (11.341) ), we arrive at the equality w(r, cj) = 
J^^ip{r],u)gk{r,r])dr]. Denoting j)(a:, cj) = p{x,u;u), we deduce from the latter state- 
ment the unique solvability of BVP ( [rm -( [rT2B . 

Now let us prove the representation l{(p) = l{(p). By virtue of (11.561) and (11.751) , 



- — ■ — - m „ 

k=i "^^^ 

m ^ „ 

= / / ^li^)9k{x,ri)p{ri)dT]yk{x]^,^k{^))dx + c 

k=l 
m „ „ 

" X] / / ^ki^)9k{r,x)yk{T](p,^k{uj))p{x)dTdx^c. (1.127) 

k=l 

The function p{x,uj;u) := p{x,uj) is a solution to BVP (11.1221) — (11.1251) with u = u, 
therefore G H^{flii \ Q) the following identity holds 



a*{p{-,uj),6{-)= / {Vp{x,io)Ve{x))-k'^p{x,co)e{x))dx 



f M^^^p{-,uj)e{-)dTR 



I Xnk{x)rl{r)gk{r,x)[yk{r,u,^k{uj)) 



(f{r],iu)gk{r,r])dr] 



dreixdx. (1.128) 
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Setting in this identity 6{x) = p{x) we obtain 



a* {'p{- , iu) , p{-)) = / {V'p{x,io)Vp{x)) — k ^{x,^) p{x)) dx 



Mf^p{-,u)p{-)dTR 



I yZ Xnk{x)rl{r)gk{r,x){yk{r,if,^k{i^)) 



/ if{r],uj)gk{r,r])dr]\ dTp{x)dx. (1.129) 



Since satisfies (11.801) — ( [1.831) and consequently (11.1051) . we have 



a{p,e)= / {Vp{x)Ve{x)-k'p{x)e{x))dx- Ml^'pedTn 
JurXh Jtr 



/ xnMQ^'z{x)e{x)dx+ Q:,he{-)dr yeeH\nR\n), 

'nR\n Jt 



which yields 



a{p,p{-,uj) = I {Vp{x)Vp{x^(jj) — k'''p{x) p{x^(jj)) dx — I p{- ^ lo) dT ji 

JnR\n Jtr 



'nR\n 



Xno{x)Qi^z{x)p{x^^) dx ^ J Q2^zp{-,u)dV. (1.130) 



Since a*{p{-]uj),p) = a{p,p{-; u)), (11.1291) and (11.1301) imply 



/ / rl{T)gk{T,x){yk{T]^,ik{^))- I (p{r],uj)gk{T,r]) dr]\ dTp{x) dx 



/ xnf,{x)Q^^p{x,u)z{x)dx + / Q2^p{-,^)zdT. (1.131) 
'nR\n Jt 



Equating (11.1271) and (11.1311) . we obtain 



l{(p)-c-^ / rl{r)gk{r,x) (p{r],u)gk{r,r]) dr] dTp{x) dx 

k=l ^ ^fc '^^k 'J ^k 



'nR\n 



Xno{x)Qi^z{x)p{x,uj) dx + J Q2^zp{-,uj) dV. (1.132) 



30 



Next, since (p{-.,cj) and z satisfy, respectively, equalities ( 11.1181) — ( 11.1211) and ( 11.761) — ( fl. 791) 
these functions satisfy also the identities 



a{(p{-,uj),6) = / xnoix) (Qi^p{x,uj) + fo{x))e{x)dx 



{Q^'pi-,cj) + go)edr (1.133) 



and 



a%z,e) = / {Xcoo{x)k 

-^Xiikix) rl{r])gk{r],x) I gk{r],<.)p{<.) d<. dr]\ e{x) dx. (1.134) 
Setting in ( 11.1331) 6{x) = z{x) and in ( 11.1341) 6{x) = Lp{x^u) and taking into notice 



that a*(2;, (p{-^ u)) = a((^(-, oj), z), we have 



Jt Jnn Jr 



XLuo{x)k{x)^{x,uj))dx 



m „ „ 

^Xnui^l) I rl{T)gk{T,x) / gk{T,T])(f{r],u)) dr]dTp{x) dx. (1.135 

J ^k J ^k 



Representation ( 11.1131) follows now from ( 11.1321) and ( 11.1351) if we take into account 

(EZSI). □ 

Remark 1. // we define a minimax estimate (p{x,uj) of the unknown solution <-p{x) 
of BVP ^.49\) - [1.5S\) as the estimate linear with respect to observations ( fi.561 ). which 
is determined from the condition of minimum of the maximal mean square error of 
the estimate taken over sets Gq and Gi, then it may he shown that, under certain 
restrictions on Gq and Gi, this minimax estimate of (p{x) coincides with the function 
(p{x,u) obtained from the solution to problem ( 11.1141) — ( 11.1211) . 
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1.5 Minimax estimation of the right-hand sides of equahties that enter 
the statement of the boundary value problem. Representations for 
minimax estimates and estimation errors 

The problem is to determine a minimax estimate of the vakie of the functional 

1{F)= / k{x)f{x) dx+ IhgdV (1.136) 
Jvlq Jr 

from observations ( 11.561) in the class of estimates linear with respect to observations 

1{F) = Y^ / Uk{x)ykix) dx + c, (1.137) 
k=i "^^fe 

where Uk G L^(ilfc), /c = l,m, c G C, and Iq G L^(Qo) and li G L'^{T) are given func- 
tions, under the assumption that F := (/, g) G Gq and the errors ^(•) = (^i(-), • • • , ^m(')) 
in observations (11.561) belong to Gi, where sets Gq and Gi are defined by (11.531) , (11.541) , 
and (11.551) , respectively. 

Definition 1.2. The estimate of the form 



- — m „ 

KF) = Y. M^yk{oc)dx^c, (1.138 



will he called the minimax estimate of 1{F) if the element u = (wi, . . . , Um) G H and 
number c G C are determined form the condition 

sup ni{F) - 1{P)\' ^ inf 

Here 

m „ 

KF) = ^ Uk{x)yk{x) dx + c, (1.139) 

k^i 

yk{x) = J^^gk{x,y)(p{y) dy + ^k{x)^ x G ilfc, k = l,m, and (p is a solution to EVP 
(11.491) — ( Tl. 5 21) at f = f and g = g. The quantity 

a := { sup E[l{F) - l{F)]^y/^ 

(/,g)eGo,|eGi 

will be called the error of the minimax estimation of expression (11.1361) . 
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Lemma 1.3. Finding the minimax estimate of 1{F) is equivalent to the problem of 
optimal control of a system described by the BVP 



z{-,u)eH\{nR\n),A), 



m 



Az{x;u) + k'^z{x;u) = ^Xnk{x) / gk{r],x)uk{r]) dr] in Qr\Q, 

k=i ^^'^ 



dz{-; u] 
du 



on r, 



dz{-,u) (2) 

— — — = Ml 'z[-] u) on Tr 



(1.140) 
(1.141) 

(1.142) 
(1.143) 



with the quality criterion 



I{u) := I Qi^{lo + z{-;u)){x){Iq{x) -\- z{x;u)))dx 

+ {Q2\h + z{--u))){h + z{--u))dT + y / rf{x)\uk{x)\^dx ^ mi . (1.144) 
Proof. Taking into account (11.1361) . (11.1371) and (11.561) . we obtain 

l{F)-l{F) = J lo{x)f{x) dx + J hgdV-^ j Uk{x)yk{x; ip, ik) dx - c 



= / lo{x)fdx+ hgdT-^ / Uk{x) / gk{x,y)(p{y) dy dx 

JQq JT k=l <^^fc 



/ Uk{x)i,k{x) dx- c 



= / lo{x)f{x) dx + hgdT-^ / gk{y,x)uk{y) dy (p{x) dx 



/ Uk{x)ik{x) dx- c 



= / lo{x)f{x)dx+ I hgdV-S^ / (xi^fe(3^) / gk{y,x)uk{y) dy\(p{x) dx 



/ Uk{x)ik{x)dx - c. 



(1.145) 
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For any fixed u = (wi, . . . , Um) G H introduce the function z{x; u) as a unique solution 
of problem ( 11.1401 ) — ( II .MSI ). According to the equivalent variational formulation of this 
problem it means that z{x] u) satisfies the integral identity 

e)= / {Vz{x]u)Ve{x)-k^z{x]u)e{x))dx- Mf\{-]u)edTR 
Jvtn\vt Jtr 

= -/ _^Xi^.(2:)/ gk{y^x)uk{y) dye{x)dx \/eeH\^R\0). (1.146) 



Set 6 = (p m ( 11.661) to obtain 



a*{z{-]u)^(p) = I {Vz{x]u)V(p{x) — k^z{x]u)(p{x)dx — / M^^ z{-]u) (p dT ^ 

P m „ 

= y^Xi^fcl^^) / gk{y,x)uk{y)dyCpdx. (1.147) 



On the other hand, since is a solution of problem ( 11.491) — ( [1.521) with f = f and 
g = setting i/j = (p and 6 = z{-; u) in ( 11.191) , we find 

afe.(.;«))= / (Vv3(,x)V;(^-*V(.x)»<i- / M<"v35M)dr« 

JnR\VL Jtr 



= / f{x)z{x;u) dx + / ^2;(-;M)(ir. (1.148) 



By Lemma 1, a*(z,(^) = a{ip,z). This identity together with ( 11.1451) , ( 11.1471) , and 
( ir™ imply 



/(F)- /(F) = / /o(x)/(a:)(ix+ / /i^dF 



^ TO /. TO „ 

/ _^Xnk{x) gk{y,x)uk{y)dyCp{x)dx I Uk{x)ik{x) dx 



= / Iq{x) f {x) dx / ligdT 
Jrio Jt 



+ / ]{x)z{x\u)dx ^ \ ^z(-; w) c/r — / Uk{x)ik{x) dx — c. 
7 rig 7r J fife 

= / f{x){lo{x) + z{x;u))dx^ / ^(/i + 2;(-; w)) (iP — / Uk{x)^k{'x) dx — c. 
J rig Jr ^ ^fe 
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The latter yields 



E 



1{F) - 1{F) 



/ f{x){lQ{x)-^z{x]u))dx-^ / g{li + z{-] u)) dr 
'fio Jt 



+E 



m „ 

/ Uk{x)ik{x)dx 



Therefore, 



inf sup E|/(F)-/(F)|2 = 



= inf sup 



/ f{x){lo{x)-^z{x]u))dx-^ / ^(/i + ^(•; w)) — c 



/ Uk{x)ik{x)dx 
k=i "^^fc 



+ sup E 

Beginning from this place, we apply the same reasoning as in the proof of Lemma 1.2 
(replacing z{x, u) by Iq{x) + z{x^ u)) to obtain 

inf sup ni{F)-l{P)\' = I{u), 

where I{u) is determined by formula ( 11.1441 ) for 

c= fo{x){lo{x) ^ z{x]u)) dx ^ / gQ{li ^ z(-]u)) dr. 
Jrio JT 

The following result follows from this lemma, 
Theorem 1.3. The minimax estimate of 1(F) has the form 



□ 



where 



c 



. m „ 

/(F) = ^ / Uk{x)yk{x) dx + c, 

= / {z{x) + lo{x))fo{x) dx + / (z + /i)5roc?r, 
jQn Jt 



1.149) 



(1.150) 



(1.151) 



and functions z and p are determined from the solution to the following problem: 

zeH\{nR\n),A), (1.152) 
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m 



Az{x) + k'^z{x) = ^Xnk{x) / gk{'n,x)uk{r]) dr] in UrX^, 



dz 

— — = on r, 
ov 

— = on Tr, 

peH\{nR\n),A), 

(A + k^)p = xnoQi\z + /o) m \ 
— = Q2^(^ + ^i). on r, 



— = Ml'p on Tr. 
ov 



1.153) 

1.154) 

1.155) 
1.156) 
1.157) 
1.158) 



;i.i59) 

Problem ( 11.1521) — ( 11.1591) is uniquely solvable. The restrictions of the solutions of this 



problem on ^^^min{i?i R2} corresponding to R = Ri and R = R2 coincide with the 
solution corresponding to R = min{i?i, i?2}- 

Estimation error a is determined by the formula a = l{Py^'^, where P = {Qi^{lo + 
z),Q^\li + -iDz)). 

Proof. Similarly to the proof of Theorem 2.1, we will show that the solution to the 
optimal control problem ( 11.1401) — ( 11.1441) can be reduced to the solution of the equation 
system ( 11.1521) — ( 11.1591) . To this end, note that there exists the unique element u G H 
at which functional ( 11.1441) attains its minimum, namely, I{u) = infy£ij/(w). In order 
to prove this statement, represent I{u) as 

I{u) = i{u) + L{u) + c, (1.160) 

where 

I{u) = / Qi^z{x]u)z{x]u) dx / Q2^z{-]u)z{-]u) dr ^'^^ / r^^{x)\uk{x)f dx 

is a quadratic functional in space H corresponding to a semi-bilinear continuous Her- 
mitian form S 



"Its continuity is shown in the course of the proof of Theorem 1.1. 
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m 



+ ^ / r,^'^{x)uk{x)vk{x)dx 
k=i '^^^ 

on H X H which satisfies the inequahty 

/ \ II 1 1 

I{u) > a\\u\\fj^ \/u ^ H, a = const, 



(1.161) 



(1.162) 



L(w) := 2Re / Qi z{x;u)Iq{x) dx + 2Re / Q2 z{-]u)lidr, 
is a Hnear continuous functional in H , and 



2= / Qi lo{x)lo{x) dx + / Q2 lihdr. 

Jilo JT 

This statement yields (see page Elj) the existence of the unique element u G H such 
that 

/(^)= inf /(^). 
Therefore, for any r E R and v E H the relations 

■^I(u + Tv) =0 and -^liu^iTv) =0 (1.163) 
dr T=o dr r=o 

hold. Taking into account that functions z{x] u + rv) and z{x; u + irv) can be written, 
respectively, as z{x; u + rv) = z{x; u)+tz{x; v) and z{x; u + irv) = z{x; u)+irz{x; v), 
where z{x; v) is the unique solution to problem (11.901) — ( Tl. 931) at u = we deduce from 
( jLTHHl ) that 

(QrH^o + z{-; w)), 5(-; v))L2(no) + (Q2 ^(^1 + ^(s ^)). ^))r 

TO 

Introduce the function p(x) £ H^({Uji \ Q), A) as the unique solution to the BVP 



(A + k^)p{x) = xno{x)QT\z{-, u) + k){x) in \ n, 
dp 



du 



= Q2 ^{z{-; u) + li) on r, 
— = Ml>p on Ti?, 



(1.165) 
(1.166) 
(1.167) 
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or to the equivalent variational problem 
a{p,e)= f {VpVe-k^pe)dx- f M^^^pOdTn 

= [ xno{^)Qi\z{-,u)+io){x)edx+ [ Q,\z{-,u)+h)edrR yeeH\nR\n). 

JnR\Ci Jtr 

(1.168) 

Setting in ( 11.1681 ) = z{-;v), we obtain 

JniiXn jvr 

= xno{x)Qi^{z{-;u) + lo){x)z{x;v)dx + Q2^{z{-;u) + li)z{-;v) dTR. 
JnR\n Jtr 

(1.169) 

Taking into account the fact that z{-;v) satisfies variation equation (11.271) with 
iIj = z{- : v) and f{x) = -J^'k^iXCiki^) J^^ gk{y, x)vk{y) which is equivalent to BVP 
(11.1401) — (11.1431) with u = V and setting in (11.271) 6 = p, we have 



a*(5(-; -u), J?) = / {'Vz{x]v)Vp{x) — k'^z{x;v) p{x)) dx — / M^'' z{-]v) pdTji 
JnR\n Jtr 



p m „ 

/ ^Xnr,{x) I gk{y,x)vk{y)dyp{x)dx. (1.170) 



Since a*{z{-] v),p) = a{p, z{-] v)), we obtain 

/ Qi^{lo^ z{-]u)){x)z{x]v)dx^ / Q2^{li + z{-]u))z{-]v) dVR 



From dUMD it follows that 



- / y'xi7fe(2^) / gk{y,x)vk{y)dyp{x)dx 
= I i gk{x,y)p{y)dy\vk{x)dx. (1.171) 



{Qi ^{lo + z{-; w)), /o + z{-; v))L2(no) + {Q2 ^{h + z{-; w)), z{-; v))t 

m p 

= -^/ r^'^{x)uk{x)vk{x)dx. (1.172) 
k=i 
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Relations ( 11.1711 ) and ( 11.1721 ) imply 

m „ 

/ rl'^{x)uk{x)vk{x)dx 
k=i ^^fe 

ni r / r 

= E 

Hence, 



Uk{x) = rl{x) I gk{x,y)p{y)dy. 
Now we can determine estimation error a. Substituting 

Uk{^) = rli^) / 9k{x,y)p{y)dy, k=l,m, 

to the formula I{u) = cr^, in which I{u) is calculated according to ( 11.1441 ) we obtain, 
taking into notice that z{t) = z{t; w), 

= Qi\k + z){x){lo{x) + z{x)) dx + / Q2\li^ z){li + z) dr 



m 

E 



9k{x,y)p{y) dy 



dx. 



Transform the sum of the first two terms. To do this, make use of equalities ( 11.1561 ) — ( 11.1591 ) 
yielding 

a{p,z)= xnoix)Qi\lo + z){x)'z{x)dx+ Q2\li + z)zdr, (1.173) 



so that 



a =a{p,z)+ I Xnn{x)Q^ {Iq + z){x)lo{x) dx ^ / Q2 {li^z)lidr 
JnR\n Jt 



+ E 



k^i ^ ^'^ 



9k{x,y)p{y) dy 



dx. 



(1.174) 



However, because z satisfies ( 11.1521 ) — ( 11.1551 ), the following integral identity holds 



m 



a*(z,e) = - _^Xnk{x) gk{y,x)uk{y)dye{x)dx ^9 e H^{Qr\Q); 
setting in this identity 6 = p, we have 



a*{z 



.P) = - y^Xr!fe(3:) / gk{y,x)uk{y)dyp{x)dx. 
JnR\n f^^^ JQk 
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The last formula, the relation a*(z,j9) = a(j9, z), and ( 11.1741 ) give us 



(7 = / k{x)Q^\k + z){x)dx / /1Q2 (^1 + ^) 



m 



/ 9k{y:x)uk{y)dyp{x)dx + ^ 

k=l k=l 

Repeating literally the end of the proof of Theorem 2.3 we obtain 



9k{x,y)p{y) dy 



dx. 



E 

k=i 

m 

E 

k^i " 



gk{y,x)uk{y) dyp{x) dx 



-2/ N 

ry I ry* 1 

ri^[Xj 



9k{x,y)p{y) dy 



and finally 



a'= / lo{x)Q^\lo^z){x)dx hQ2\h^z)dr = l{P). 
Jfio Jr 



□ 



In Theorem II .41 stated below we obtain another representation for minimax estimate 
/(F), not depending on the form of functional /. 



Theorem 1.4. The minimax estimate of 1{F) has the form 

1{F) = /(F), 



(1.175) 



where F = (/, g), f{x) = f{x, lj) = ^p{x, u) + fo{x), g = g{-,uj) = ^joPi': cu) + 
go, and function p = p{-,u) is determined from the solution to problem (11.1141) — (11.1211) . 



Proof. Taking into notice ( 0491 ) - ( ITml ) . we find 

/(F) = ^ / Uk{x)yk{x) dx + c 
k^i ^^'^ 

= X] / / rl{T)9k{r,x)yk{T)p{x)dTdx + c, 
where c is determined from (11.1511) . 



(1.176) 
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Next, repeating literally the proof of Theorem 1.2 on page EO, we arrive at the 
relatioship 



a* (p{- , uj) , p{-)) = / {'Vp{x,uj)'Vp{x)) — k'^p{x,cu) p{x)) dx 



Mf^p{-,u)p{-)dTR 



/ / X^X^yl{'r)9k{T,x)[yk{T) 

(p{r], u)gk{r, r])dr]^drp{x) dx. (1.177) 
Taking into account that p{x) satisfies (11.1561 ) — ( 11.1591) . we obtain 



a{p,e)= / {Vp{x)Ve{x)-k^p{x)e{x))dx- Ml^^pOdTR 
JnR\n Jtr 



= / xno{x)Qi{z + io){x)e{x)dx+ Q,'{z + h)e{-)dr ye e H'{nR\ny, 

JnR\n Jr 
consequently, 

a{p,p{-jU)) = / {Vp{x)Vp{x^Lo) — k^p{x) p{x^uj)) dx — / M^^}^ p p{- ^ uj) dV ji 
JnR\n Jtr 

= f x^,{x)Qi\z + k){x)p{x,uj)dx+ f Q^\z^h)p{;Lu)dr. (1.178) 
JnR\n Jt 

From (11.1771) and the latter equality, it follows that 

XI/ / ^ki^)9k{r:x) lyk{r) - ip{r],Lj)gk{T,r]) dr]j dTp{x) dx 

= / Xn,{x)Q^^p{x, u) {z{x) + /o(x)) dx + / Q^^p{-,uj){z + h) dV . (1.179) 
JnR\n Jr 

Equating f^JTE) and f fTTTm , we find 

. -~- — ^ Tn n n ft 

l{F)-c-^ / rl{r)gk{T,x) (p{r],u)gk{r,r]) dr]dTp{x) dx 

k=l "J ^k *J ^k 

+ / xnS^)Qi\{x)p{x,uj) dx ^ I Q2^lip{-,u) dV 
JnR\n Jr 
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= / xno{^)Qi^z{x)p{x,u) dx ^ / Q2^zp{-,u) dT. (1.180) 
JnR\n Jt 

Next, since <-p{-^co) and z satisfy, respectively, equalities ( 11.1181 ) — ( 11.1211 ) and ( 11.1521 ) — ( 11.15' 

these functions also satisfy the identities 



a{(p{-,u),9) = / xno{x) {Qi^p{x,oj) + fo{x))9{x)dx+ / (Q^'^pi-,^) + go) 9 dT 
JriiiXfi Jr 

(1.181) 

and 

p Tin n n 

a\z,9) = -j _^Xnk{x) rl{r])gk{'n,x) gk{r],<;)p{<;) d<; dT]9{x) dx. (1.182) 

Setting in ( 11.1811 ) 9{x) = z{x) and in ( 11.1821 ) 9{x) = (p{x^uj) and taking into account 
that a*(2;, (p{-^ lo)) = a((^(-, cj), z), we obtain 

JT Jno Jr 



'nR\n 

m 



IIV n n 

^XnAv) i rl{T)gk{T,x) I gk{T,r])(p{r],uj))dT]dTp{x)dx. (1.183) 



From ( fTTM ) and ( ITTK^ . it follows, in view of ( fTTni ). that 



1(F)— / z{x)fo{x)dx— / zf/oo^r- / Iq{x) fo{x) dx — / /i^^oG^r 

+ / 2;(a:)/o(a:) (ix + / zgodT 
Jfio Jt 

= / /o(x)Q;fV(2^? ^) '^a: + / liQ2^p{-,io) dV, 
Jrio Jt 



thus 



1{F)= / /o(x)(/o(x) + Qr'l^(^,^))^^+ li{h + Q2^p{-.^))dT 
Jno Jt 



= / /o(a:)/(x)(ix+ ligdT = l{F). 
Jnn Jt 



□ 



Remark 2. // we define a minimax estimate F(x,uj) of the element F = {f,g) 
as an estimate linear with respect to observations (il.56[) . which is determined from 



^Here / and g are the functions entering the statement of BVP (|1.49p - (|1.52p and F = (/, g) € Go- 
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the condition of minimum of the maximal mean square error of the estimate taken 
over sets Gq and Gi, then it may be established that, under certain restrictions on 
Gq and Gi, this minimax estimate of F coincides with the element F = (f^g), where 
f = Qi^p{x,uj) + /o(a:) and g = Q2^'Jdp{-,^) + go, and the function p = p{-jUj) is 
determined from the solution to problem ( 11.1141 ) — ( II ■1211 ). 

Using Theorems 1.1 — 1.4 together with the solution techniques employing the so- 
called DtN finite-element methods elaborated for problems ( 11.141) - ( [1.161) and ( 11.411) — ( [1.431) . 
one can construct algorithms of numerical solution to problems ( 11.761) — ( [1.831) . ( 11.1141) — ( 11.1211) 
and ( 11.1521) — ( 11.1591) and obtain the required minimax estimates. 

Remark 3. All results of this section remain valid in the three-dimensional case. For 
example, finding minimax estimates of the solutions to the EVPs for the Helmholtz 
equation that describe diffraction of acoustic waves by an obstacle Q G can be reduced 
to the solution of integro- differential equation systems ( 11.761) — ( 11.831) and ( 11.761) — ( 11. 831) : 
the domain \ Q, should be replaced by \ Q, and plane domains Qj, i = 1, /c, on 
which observations are made should be considered, as well as supports Qq and ujq of 
functions f and /q, as spatial domains. The Sommerfeld condition 

- ikif = o(l/r^/^), r =\x\ = + r ^ oo (1.184) 

should be replaced by 
dip 



dr 

defined on a sphere Tji of radius R 



ikip = o(l/r), r = |a:| = x\ + x^^ Xg, r oo; (1.185) 



operators Mj^ il) defined on a circle Tji should be replaced by the following operators 



{Mi^^^^J){R, 0) := ^ E -^Wr^ E ^rnnYUO. 0), (1-186) 

^ 71=0 [kR) n=-r 



where (i?, 6, </)) are the spherical coordinates of the point x = (xi, X2: x^) G F^^, Umn = 
Jq Jq^ tlj{R, 9', (j)')Ymri{(^',(j)' ) sin 9' d 6' del)', hn\x) are the spherical Hankel functions 
{j = 1,2), Y^n{9^(j)) := l^f^^^^^^Pn"'\cos9)e'^'^ are the normalized spher- 
ical functions, and P^{t) are the associated Legendre functions (—n < m < 
n = 0, 1, . . . , ooj; and formula ( 11.181) should be replaced by 

(/p(rp, 9p, 0p) = ^ ^ TbrrSr E ^-^nY^n{9p, rp > R. (1.187) 



(kR) m=-T 
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Remark 4. The analysis performed in Section 2 enables us to state that all the results 
of this section remain valid, for example, when in BVP ( 11.491 ) — ( [1.521 ) the Helmholtz 
equation ( 11.501 ) is replaced by 

-{A + Pn{x))if{x) = f{x) in \ Q, (1.188) 

where n = 2, 3, A; = const > 0, the function n G C{W^ \ Q) is positive in M"" \ Q, and 
n{x) = 1 in the domain ^Ir^ \ for a certain Rq > 0. 

In this case one has to choose in the equation systems ( 11.761 ) — ( 11.831 ) and ( 11.1141 ) — ( 11.1211 ) 
which specify minimax mean-square estimates the value of R greater than Rq, setk = 
and replace in equations ( 11.771 ). ( 11.811 ) . ( 11.1151) . and ( 11.1191) by k'^n{x). 

Note also that applying DtN finite-element methods to problems ( 11.761) — ( [1.831) and 
( 11.1141 )- ( 11.1211) one can construct approximate methods of their solution. 
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PART 2 



Minimax estimation of the solutions to the boundary value prob- 
lems from observations distributed on a system of surfaces. Reduc- 
tion to a surface integral equation systems 

2.1 Statement of the problem 

Before to formulate the estimation problem which is the subject of analysis of the 
present chapter, let us introduce the necessary notations and functional spaces 

Let A be closed or unclosed (n — l)-dimensional Lipschitz surface in R*^. By dA 
we will denote the element of measure on surface A and by i^^(A) the space of square 
integrable functions on A. 

Let 7 be an unclosed (n — l)-dimensional smooth C°°-surface in its boundary 

whose points do not belong to 7, d^Ci^ = 0, and 7 a closed smooth C°°-surface ((n — 1)- 
dimensional manifold without edge) that contains 7, 7 C 7, and divides M.^ into two 
domains, bounded and unbounded. Set 

^1/2(7) := {^1^ : V e H^'\l)]. 

The norm in space i7^/^(7) is determined according to the formula 

Denote by -ff "''^(t) = {-ff ''^{7))' the space dual to H^^'^{'y). Below, the duality relation 
< r^w >j on i7^^/^(7) x i7^/^(7) will be also denoted by J^rwd-j because for this 
relation the condition (*) on page M is valid in which F should be replaced by 7. Note 
that the elements of i7~^/^(7) extended to 7 \ 7 by zero values belong to H^^^'^{'j). 

Let p{x) be a function regular on surface 7 which is equivalent to the distance from 
a point X to the boundary of 7 (this distance will be denoted by d{x^d^)) in a 
vicinityj of d^. Following [30] and [SO] introduce the space 

Hll,\^) := {u e H'/\j), p-''\ G L2(7)} = {u^ H'/\^), u G H"\^)} 



®It means that Ximx^xo dix'^d-y) ^ ^ const ^ Vxq G d"f. Such functions exist because is an infinitely 
differentiable manifold 1301. 
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with the norm 

M M (w m2 m -1/2 ii2 

\\'A\ull,\,)=y}\^\\ll-l-{,)^\\p ^lli^(7)J ' 

where function u denotes the extension of w by on 7 \ 7. 

^00 (7) ) we will denote a space conjugate to Hqq (7). Then, in line with [5^, 
p. 43, we have 

(^o^o'(7))' = {/ = /o + /i, /o G ^-'/^(t), P'^'fi e i^'(7)}. 

Now let us formulate the estimation problem. Assume that the state <-p{x) of a 
system is determined as a solution to the Neumann problem^ 

^eHl^^{{R'\n),A), (2.1) 

(A + = in R^\n, (2.2) 

^ = h on r, (2.3) 

dip 

— ikif = o(l/r), r = |a:|, r 00, (2.4) 

whercl G is a bounded domain with a connected complement such that 80, = F 
is a surface of class C'^, h G ^^(r). 

It is known that problem (12.11) — ( 12^ is uniquely solvable. 

Let 7j, = 1, A/^, be smooth simply-connected oriented surfaces in with smooth 
boundaries d^i, djiO^i = 0, contained in the domain \ H that have no intersections 
pairwise, 7^ fl 7^ = 0, i 7^ j, 7^ C \ Q. Let the orientation of 7^ be determined by a 
continuous family of normals ^{x), x G 7^. 

Assume that on surfaces 7^ the following functions are observed 

yP{^)= f Kl'''\x,yMy)dj,^+ [ Kl'''\x,y)^dj,^ + rjl'\x), (2.5) 
vf\^) = f Ki'^'\x,yMy)d^,^^+ f Kf^^{x,yf-^d^,^-rvf\x), (2.6) 



this chapter we will restrict ourselves to the case n = 3. The results obtained for n = ?> remain valid for n — 2 
after corresponding replacement of Sommerfeld radiation condition and fundamental solution. 

*It is known (see |40], page 221) that there exists a uniquely determined continuous operator which we denote by 
-S- and which maps space i7^(ri, A) into space 77-i/2(r) and is such that Vu £ H^{n,l\), Vw S H^{n) the following 
representation (Green's formula) holds: X]"=i/n^^^^ = ~ /n '^^ ^^2; + /pl^udF, where the integrals over T 
should be understood as the duality relations < ^,fJ-v > on i/~^/^(r) x H^/'^(r). This operator is called the normal 



du 
dxj 

the unit normal vector of F external with respect to domain fl and cos(i/, Xi) is the zth directional cosine of 



derivative in relation to —A; the operator ^ is defined by ^ = cos(z^, x^) when u € C°°(r2), where i/ is 
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where '^{x) is the solutioij^ to BVP (O) - (j2j) , r]f\x) and r]f\x) are the observa- 
tion errors that are choice functions of random fields defined on surfaces 7^; K^"^^ G 
L'^ili X 7i), T^j = 1,2, are functions defined on 7^ x 7^; and integral operators G^p 
with kernels Kp'^\^,x) defined according to 

Gl'Uix) = I Kp'\^, x)md7^,: J = 1, 2, (2.7) 
are Hnear bounded operators acting from L {'ji) to Hqq (7^), i = 1,N (as an example 
of such kernels, one may take degenerated kernels K^-'''^\^,x) = Yll^i ^r'^\0^r''\x) , 
where ai'^^ G L\j,), bi'^^ G H^i^j,)). 



From the physical viewpoint, observations of the form (12.51) , (12.61) enable one, e.g. 
in stationary problems of hydro acoustics, to observe independently both the pressure 
and the normal velocity component as well as their linear combinations on a system of 



surfaces 7^, i = 1, A^. 

Denote by Go the set of functions h E ^^^(r) that satisfy the condition 

J Ih-hofqjdr < 1, (2.8) 

where ho G L'^{T') is a given function. By Gi we denote the set of random vector- 
functions fj{-) = (r/l^^ (•),..., ?7^^(-), ?7p^(-), . . . , ?7^^(-)); their components fjl^\x) and 
ffl (x) are random fields defined on surfaces 7^, i = 1,N having square integrable 
second moments and satisfying the conditions 

Er/f ^(x) = 0, Bfif\x) = 0, 1 = TJV, (2.9) 

AT 2 ^ f 2 

rf7. + E / E|^f^(^)l' {rf\^)) dl. < 1, (2.10) 

where qi{x) , rp\x) , rf'\x) , i = l^N, are functions continuous on T and 7^, respec- 
tively, that do not vanish on these sets. 

Assume also that in equalities (12.11) — ( 12^31 ) function h{x) and the second moments 
E|?7j^^''(a:)|^ and E|?7p^(a:)p of random fields r]^^\x) and rif\x) in observations (12.51) and 
(12.61) are not known exactly, and it is known only that 

h G Go, r/(.) = (r/f ), . . . , r/^)(.), r/f )(•), • • • , V^n\-)) ^ G,. (2.11) 



^Note that for any subdomain to such that to G fl, the solution ip to problem l|2.ip — (|2.3p belongs to H'^{uj); therefore 

M I 

<i/ 1 7. 



according to the trace theorem, fl-y^, ^-^^ G L^{'-fi), i = 1,-/V, and integrals l|2.5p — (|2.6p make sense. 
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Suppose that a function Iq G L'^{coo) is defined in a domain wq, cjq C \ The 
problem is as follows: given observations ( 12 .SI ). ( 12.61 ) of the state <-p{x) of a system 
described by the Neumann BVP ( 12.11) — ( |2r3l ) under the conditions that h E Gq and 
?7 G Gi, estimate the value of the linear functional 

l{ip)= / lo{x)(p{x) dx (2.12) 
in the class of estimates linear with respect to observations that have the form 

N „ 

^) = T. U\^)yl'\^) + uf\x)yi'\x)) + c, (2.13) 

where \ uf^ G ^^(7,), i = Tj^^ c G C. 

Put u := {u^^\ u^i\ uf\ u%^) eH:= {L\^,) x . . . x L\^^)f . 



Definition 2.1. An estimate 

^) = Y. [^P{x)yf\x) + uf\x)yf\x)) d^, + c, 

i=l 

in which functions uf'\ uf^ and number c are determined from the condition 

sup E\l{^) - ^ inf (2.14) 

l^r^ -^n ueH,c£iL 



where 



N „ 

^) = J2 U\^)yl'\^) + uf\x)yP{x)) d^, + c, (2.15) 

yl'\^)= f Kl'^'\x,yMy)dl^,+ I Kl'^'\x,y)^d^,^ + f)l'\x), (2.16) 

yf\^) = I Kf^'\x,y)^{y)d^,^+ I Ki'-'\a:,y)^dj,^ + fjf\x), (2.17) 



a: G 7j, i = l,N, 

and (p{x) is the solution to the Neumann BVP at h = will be called a minimax 
estimate of expression ({2. IB) . 
The quantity 

a:={ sup E|/((^)-/(^)|2}i/2 (2.18) 
will be called the error of the minimax estimation ofl(ip). 
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2.2 Auxiliary statements 

In this section we will prove that finding the minimax estimate is equivalent to a certain 
problem of optimal control of a system described by elliptic equations with conjugation 



conditions on surfaces 7^, i = 1, A^. 

In order to state the conjugation problems under study and prove the existence of 
their solutions it is necessary to introduce the corresponding Sobolev spaces and trace 
theorems for surfaces with edges. First, let us formulate several definitions. 

For / G V{W^) the function u{x) := Jj^„ ^k{x, y)f{y)dy solves —Au — k'^u = / (and 
complies with Sommerfeld radiation conditions). Here 

^M'\k\x-y\) in M\ 



in R3 



47r \x—y\ 

is the fundamental solution to the Helmholtz operator. 

Introduce the Newton potential operator {Nkf){x) := J^ri^k{x,y)f{y)dy- 

Lemma, (see, for example, J^) Nk can be extended to a bounded operator Nk '■ 

For i/j^x^ C'(r) introduce the functions 

(VrV)(x) := j^^k{x,y)ilj{y)dTy, (2.19) 

{wh){x) := f ^^^^x{y)dry, x^r, (2.20) 

called the single and double layer potentials. Let us formulate the results contained in 
|50| in the following form. 

Lemma. Vy andWy can be extended to bounded operators Vy '■ H^^/'^CT) H^^^{M.^)n 
Hl^ (n U (M^ \ n), A) and : H^/\r) ^ H^^ (Q U (M^ \ Q), A) |^ 

Let 7o be a smooth bounded oriented simply-connected unclosed smooth surface in 
with a smooth boundary ^70, ^70 Pi 70 = 0. Let its orientation be specified by a 
continuous family of unit normals ^{x), x G 70. Denote by such a bounded open 
set in with a smooth boundary OUq To containing surface 70 that the normal 



^° Here, the boundedness of, e.g., the operator Vp : H ^^^{T) H^^^{M.^) means that for any cutoff function 
a e Cg^mp(I*^) the operator aV^ : H-'^/'^{T) H^{M.^) is bounded. 
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vector V to 70 is directed outside Hq; by 70- we denote the side of 70, whose orientation 
coincides with that of the external side of surface dQ.Q and by 70+ the opposite side of 

Let 71, . . . , 7iv be the surfaces of the type introduced above (see p. HG]) with the sides 
7i+, 7i_, . . . , 7iv+, 7iv- respectively. 

Set Q' = (R^ \ ^) \ ^f=ili si-nd for any function v defined in Q.' denote by v\^^^ and 



v\ 



the restriction of v to 7^+ and to 7^-, i = 1, A^, respectively. 



,(1) Ji) 



Assume that the functions 5^, ^ G H^q^^{^') =: L^Qj^p(il'), a, (5 G H^^/'^{r), uf\ g^^' G 
^00^(7*)' ^ G iif^^/^(7i), = 1, A/" are defined in domain il'. 

Consider two problems. 

1. Find function u that satisfies the conditions 



u 



{A + k^)u{x) = g{x) in Q', 
du 

du 

(1) 



= a on r, 



to- 



on 7 



1,A^, 



du 



d 



V 



(2) 

= tj\ on 7„ 



z = 1,A^. 



— iku = o(l/r), r = Ixl, r ^ 00, if Im/c > 0. 

or 



2.21) 

2.22) 
2.23) 
2.24) 
2.25) 



2.26) 



du I 

1 



Here = u\^^^ - Ai^^ ^ ^oo^(7*), ^It.+ '^It.- ^ H^''^{li). [%]^^ ^ 

e ^-^/'(7.), S;|,.,,S;|,._ e (//of (7.))', ^ = T:A^, tir e ^"^/^(F), and 
(12.221) — ( [2.251) should be understood as equalities of elements from spaces, respectively, 
L\^'), H-^I\T), //o'o'(7.), and H-^l\^,), i = iJJ . 
2. Find function v that satisfies the conditions 



^G//f^^(^^^A), 

- {/\ + W)v{x) =g{x) in Vl\ 
— = a on i , 

W\l^ = d'P on 7„ z = l,A^, 



(2.27) 

(2.28) 
(2.29) 
(2.30) 
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dv 
du 



Q 



(2) 



on 7^, 



dv 
dr 



+ iku = o{l/r), r = \x\, r 



(X), 



(2.31) 
(2.32) 



where [v],, = ^|,., -^|,._ G hI',\^,), v\,^^M,.- ^ H'^'M. [f],, = t 

^"'^'(7.), tl,^, , e (//o'/'(7.))', ^ =T:A^, tir e H-'/\r), and (j228])-(l23I]) 

should be understood as equahties of elements from the corresponding spaces. 

Remark 5. The choice of spaces for problems 1 and 2 is governed by the trace theorems 
(see lEU^, pp. 44, 45 and gU^, pp. 180, 181). 

In order to prove the existence and uniqueness of solutions to problems 1 and 2 we 
formulate one more known result. Namely, let 70 and 70 be the surfaces introduced on 
p. EOl By virtue of the definition of spaces iifg 0^(70) and i7~^/^(7o), the elements of 
these spaces extended by zero on 70 \ 70 are the elements of iif^/^(7o) and H^^/'^{^q), 
respectively. 

Set 

Xl := {u G L»'(M^), u G H\nR \ 70), Au G H\nR \ 70) V sufficiently large R}. 



dv_\ ^ 



(2.33) 

If the tilde sign marks the zero extension on 70 \ 70 of an element defined on 70, the 
following statement holds. 

Theorem 2.1. Let p G //oo^(7o) and p' G H-^/'^{-fo). ThenW^^p := W^^p andV^^p' := 
V^^p' belong to X^^, and V p' G E 



70' 



Formulate the BVP: find u G i7j^Q^(IR^ \ 70) satisfying 



M70 = 



du 
dr 



R'\7o), 




(2.34) 


= 


in IE 


\ 70, 


(2.35) 


du 
du 


70 


P, 


(2.36) 


r - 


= \x\ 


, r 00. 


(2.37) 



The next statement is an immediate corollary of the last theorem. 
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. Then BVP (2J^-(2M) has a 



1 /o 

Theorem 2.2. Let p G Hqq (70) and p' G H~ 
unique solution u G X}^^ which, for x outside 70, can he given by the formula 



(2.38) 



Here [u]^^ = u 
G H 



7o^ 



u\^^_ G HJq (70), wUn, , u 



du 
dvA 



70- 



du 
dvA 



70+ 



dvA 



70- 



7o+> -no- 6 H^'\l<3), 
6 f^^o^'(7o) 



du 




du 


_dvA_ 


70 


dvA 



70+ 



Now let us prove, e.g. for problem 2, that the following statement is valid. 
Theorem 2.3. BVP ^2.27\) — f2.3^) is uniquely solvable and the estimate 

\\v\\m{n'nnR) < C'o(||^IU-i(r!o) + ll«llff-i/2(r) 



N 



N 



* 1 1 ^00' (7.) 



+ E 11^?^ 11^-/^(7.))' (2-39) 



holds, where Qq is the support of function g and Cq is a constant which does not depend 



on g, a, g^^^\ and g)^\ i = 1, N 



(2) , _ 



Proof. Set 



N 



ViiX 



N 



i=l 



where, in accordance with definition on page EH and W"'^ are single and double 
layer potentials determined on unclosed surfaces 7^, ? = 1, A^, corresponding to the wave 
number —k. Then by theorem 12.21 function vi{x) is the unique solution to the problem 



- {A + k'^)vi{x) = in W, 
— = on r, 

OP 

(2) 



on 7,, i = l,N, 



du 



Q 



on 7 



11 



hN. 



dvi 
dr 



+ ikvi = o(l/r), r = |x|, r 00, 



2.40) 

2.41) 
2.42) 
2.43) 
2.44) 

2.45) 
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Here [v,],„ = - t.i|„_ e ff„'/'(7,), e ^'/^(t.), [f^],, = f^l 

Denote by the unique solution to the problem 



(A + P)z;2(x) = 0, 



du du 

dv2 



on r, 



dr 

and by the unique solution to the problem 



+ ikv2 = o(l/r), r = |x|, r ^ oo, 



2.46) 

2.47) 
2.48) 

2.49) 

2.50) 

2.51) 
2.52) 

2.53) 
2.54) 

V^J'gf\x) and W~J^gf\x) = W^^^f ^(x), and boundedness of operstors V^^ and 
in the corresponding spaces, we obtain the following estimates for vi 



dp 

dv3 



= a on r, 



dr 



Then the function 



+ ikv^ = o(l/r), r = \x\, r oo. 



v{x) = vi{x) + V2{x) + V3{x) 



will be the unique solution to problem ( 12.281) — ( [2.321) . From equalities V^/^^?,*"^^ 
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N N 

IL, II . ^ ll;^(2)| 



lbi||Hi(f^'nf^fl,A) - C'l^ll^/ ||^i/2(^^), |bi||Fi(i7'ni7fl,A) - ^2 ||g. 

i=l i=l 

and hence 



^1 



lki(f^'nf^fl,A) ^ ^3 I X] Ik/ lli/V2(^^) + 2^ 11^/ ) • (2.55) 

\i=l i=l 

For functions V2 and -^a we apply estimate ( 11.341) to obtain 



MlminnXn) < C'4||7Af^i||ff-i/2(r), (2.56) 



'^Here and below C'i are constants that do not depend on the data of the problems in question. 



53 



< C's {\\9\\H-^{ni^\n) + Il«lli7-i/2(r)) • (2.57) 
Using the trace theorem ([IQ], pp. 180, 181) and ( 12.551 ). we prove the estimates 



N 



N 



(i)i 



i=l 



(2.58) 



From (|234l)-(|238D it follows that 

Iklliji(i7'nf}/j) < Co(^\\g\\H-i{n'nnR) + ||<5||7j-i/2(r) 



(2) I 



i=l 



i=l 



□ 



For problem 1 the corresponding theorem is obtained in a similar manner. 



2.3 General form of the guaranteed estimates and expression for the esti- 
mation error 



Introduce, for every fixed u ^ H = (L^(7i) x . . . x L^(77v)) the function z{x; 
solution to the problem 

- {A + P)z{x;u) = Xuja{x)lo{x) in Q\ 
dz 



= on F, 



u j as a 

(2.59) 
(2.60) 
(2.61) 



[z{x;u%^ = 



dz{x; u) 




dvA* 


7, Jl^ 



dz{x] u) 
dr 



on 7,, i=l,N. 
+ ikz{x; u) = o(l/r), r = |a:|, r oo. 



(2.62) 



(2.63) 



Lemma 2.1. Finding the minimax estimate of l{(p) is equivalent to the problem of 
optimal control of the system described by BVP Ii2.59\) — ^2. 6^) with the cost function 



I{u^i\ . . . , wS^\ u'i\ . . . , wj^^) = J Qi ^{x)z'^{x] u) dr 
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N 



+E 



r- [x 



-.1 ^7. 



N 

))-\4\x)fdj, + J2 / {rf\x))-\uf\x)fdj, ^ min. (2.64) 



Proof. Denote by Qi an open subdomain in R^\Q (Q^ C R^\Q) such that dQi contains 
7j, its boundary dfli is simply-connected and smooth, and the normal vector u to 7^ is 



directed outside Qi. We also assume that fl il^ = for all i 7^ j, j = 1, A^. 

Set := {Qn \ ^) \ U^^Q^ (i? is assumed to be sufficiently large), 7^ = dQi and 
denote by 7i_ and 7^+ the external and internal sides of surface 7^. Next, simplifying 



the notation in the surface integrals, denote by z+, ( 5^ ) and 



z 



dz 



lk+ ' dvA* 



or z 



7fe+ 



(?2 

7fc+ ' 



and 2; I 



dz 



7fe+ 



7fe- 



dz 



Ik- 



the traces 
of functions 



we 



Ik- ' dvA* 

z{x; u) or on sides 7^+ (or 7fc+) and 7fc_ (or 7^-) of surface 7^ (or 7^) 

Taking into consideration relationships (12.151) — ( [2. 171) . (12.591) — ( [2.631) . and applying 
to (p{x) and z{x]u) in domains il^, = 1,A^, and the second Green formulai^^ 
obtain, using the equalities [z{-] u)]a^.\^. = = 0, 

/(^) - m = 

^ N „ 

= / lQ{x)(f{x) dx-^ lul^\x)yP{x) + uf \x)yP{x)] d-f^ - c 



N 



N 



[x)lo{x)(f{x)dx^'^^ / XLUoi^)^oi^)'fi^) 

1=1 "^^^ 



^•^1 -'7ix 



/ (f{x){A + k'^)z{x]u)dx — '^^ / (^(a:)(A + A;2)z(a:; w) 



1=1 ^7^ 

TV 

^/ 

z=l ^7* 



(1), 

u] ix] 



(1,2), 



It 



(2)/ ^ 



(2,2), 



It 



dm 

dvA* 

dm 

dvA* 



d% 



d% 



^^One can apply the second Green formula because z(-; u), ip G H^{riii), z(-; u), ip G H^[VLi), i = 1, N, and Az, Aip e 
L'^ifln), Az, Aip e L^{fl,), i = T^/V, 
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N „ N „ 

i=l -^^i i=l -^Ti 



c 



= - / {A + k'^)(f{x)z{x;u)dx + [ - [ {A + k'^)(f{x)z{x;u) dx 



N 



N 



i=l -"^i 

N 



dv 



dv 



i=l Jli -^li 



H 



N 



-E 



dm 



(2,2). ^,(2), 



N N 

i=i "^ii 1 



— c 



=1 -^Ti 



= / z/i(ir + ^^ / u) — w)) 



T 

iV 



d(p{x] 
dv 



dz{x, u) 



du 



+ 



'dz{x, u) 
dv 



^{x) dji 



N 



/ m / Kl'^'\C.^h^\0 + Kl'''\C.^)uf\0 dj,^dj,^ + J:R{z{-;u),^) 

Af '(C) + iff ''>(5, x)«f>({)| d7i, d7i. 



-E 



d(p{x) 



. . ^-7 7-7 



- ^ Nr. ^ 



— c 



= / hzdr + ^R{z{-;u),<f) 



N „ N „ 

^ / uf\x)fjl'\x) d^,-Y ^f\^W'\^) dl^ - c, (2.65) 
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where by I1r{z{-; w), ip) we denote 



— - d(f{x) ( dz[x;u) 

^z{x;u)—^ (p{x)\^ — — — ) I dTu 



Since z{-; u) and (p satisfy, respectively, the Sommerfeld radiation conditions ( 12.631 ) and 
( 12.41 ) we obtain an estimate for Yiji{z{-] u),ip), 



T,r{z{-]u),(p) = / z{x]u] 



d(p{x] 
dR 



ik(p{x) ) dTji 



= / 0{l/R)o{l/R)dTR- f 0{l/R)o{l/R)dTR = o{l) npn R^oo. 
From here, passing to the Hmit as ^ cxd in ( 12.651 ), we obtain 



1{lp) — l{(p) = i hz{-]u)dT 



N 
z=l 



T 

N 



)fjl^\x) d-f^-Y^ / uf\x)fjP{x) d-f^ - c, 
1=1 -^y^ 



The latter equalities together with conditions ( 12.91 ) — ( 127701 ) and the known relation 
= E^^ — |E^|^ that couples dispersion of random variable ^ and its expectation 
E^, yield ^ 

inf sup E|/((^)-/(^)|2 = 



= inf sup 




hz{-] u)dT — c 



N „ N „ 

sup E ^ / uf\x)fif\x) d^i + Y^ I uf\x)fif\x)d^, 



2.66) 



In order to calculate the supremum in the right-hand side of ( 12.661 ) make use of the 
Cauchy— Bunyakovsky inequality. Introducing the notation 

y = J (h- ho)z{-]u) dr, 

we prove, using relation ( 12.81 ), the inequality 



\y\ < \ l^qr{^Mx-.u)\ dry x [J^\h 



hol'^qidT 
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in which the equahty holds at /i G Gq and 
Therefore, 



a 



(2.67) 



inf sup 



ceC 



heGo 



hz{-] u) dT — c 



at 



Similarly, 



c = y z{-] u)hQ dr. 



(2.68) 



sup E 



N 

i=l -^^^ 



)f}f\x)d-i, + ^ I uf\x)fif\x) d-i, 
■ 1 J 



1=1 

N 



< 



N p ^ r 

J2 / {rl'\x)r'\4\x)\'dj, + Yl / {r?\x))-'\uf\x)\'dj,. 



It is easy to see that in this inequality, the equality holds when f/ is a random vector- 
function with the components 

i(rf\x))-h,f\x) 



.(2)^^U-2l„,(2), 



N 



,(2)/^\\-2U,(2), 



7» 



1/2' 



1/2' 



(2.69) 



where is a random value such that = and E^^ = 1. The latter facts yield 



inf sup E|/(^) - /(^)p = I{u^^\ . . . , w^;, <^ ...,u 



,(1) .,(2) 



777-2 



where functional / is determined according to (12.641) and the infimum with respect to 
c is attained at c = z{-;u)hodr. The lemma is proved. □ 
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Solving the optimal control problem ( 12.591) — ( [2.641) . we arrive at the following state- 
ment. 

Theorem 2.4. The minimax estimate of the value of functional l{(p) has the form 

N „ 

= {4\^)y^'\^) + uf\x)yf\x)) d^, + c, (2.70) 

i=l 



where 



c = y zKq <ir, 



v!l\x) = {rY>{x)Y 
uf\x) = irf\x)y 



(1), 



MO 

du 

MO 

du 



(2.71) 



(2.72) 



and function p{x) is determined from the solution to the problem 



(A + P)z{x) = XLUa{x)lo{x) in Q', 



dz 
dv 



= on r, 



z[x)\^^ = 



dz{x 



d-fi^, 



du 



dz{x) 



-(2,1), 



dr 



+ ikz{x) = o(l/r), r = r oo, 

(A + k^)p{x) = in Q', 
dp _2 



— = z on r, 



du 



0, z = l,iV. 



(2.73) 

(2.74) 
(2.75) 

(2.76) 



(2.77) 

(2.78) 
(2.79) 
(2.80) 
(2.81) 

(2.82) 
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dp{x) 
dr 



ikp{x) = o(l/r), r = |x|, r oo. 



(2.83) 



where in ^2. 72) p{y) and denote the values of the traces of functions p and 

boundary values of its conormal derivatives on different sides of surface -ji. Also, 
^ ^^(7*)' ^ = MV- Problem fKlE)-fIM) 'is uniquely solvable. 
The error a of the minimax estimation of l{ip) is given by the formula 



a 



1/2 



Iq{x)p{x) dx 



(2.84) 



'Wo 



Note that if we replace in (12.761 ) and (12.771) ul\x) and uf\x)^ i = 1,A^, by their 
expressions in the right-hand sides of (12.721) . then these functions may be excluded from 
the equality system (ETSl) - (E^SD . 

Proof. Let us show that I{u) is a quadratic function on H. Indeed, since the solution 
z{x] u) to problem (12.591) — ( [2.621) can be represented as z{x] u) = z{x] u) +2;o(a:), where 
z{x;u) is the solution to this problem at lo{x) = and zo{x) is the solution to the 
problem 

zoeHU^',A), 

- (A + k'^)zo{x) = Xu;o{x)lo{x) in H', 
dzo 



du 



on r, 



(2.85) 
(2.86) 
(2.87) 



[^0(3^7 J 7. 

dzQ{x 
du 



on 7i, i = l,N, 



= on 7„ 1 = 1, N, 



dzQ{x) 
dr 



+ ikzo{x) = o(l/r), r = r 00, 



(2.88) 
(2.89) 



functional I{u) can be represented as 



where 



I{u) = I{u)^L{u) + J q^'\zo\'dr, 



i{u) = J qi'^{x)\z'^{x;u)\ dr 

i=l ^7? ^ 1 



=1 ^7^ 
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From inequality ( 12.391 ) and our assumptions concerning operators of the form ( 12.71 ) . 
we deduce, taking into account that || • ||^-i/2(^^) < c|| • ||l2(7,) (c = const > 0), the 
inequahty 

^ " 2 xl/2 



i=i -^^^ 



2.90) 



where ci = const that does not depend on u. Taking into account ( 12.901 ) and the 
trace theorem from [1], we see that u ^dz{-]u) is a bounded hnear operator that 
maps Hilbert space H in H^/'^^iT). From the latter statement, it follows that I{u) is a 
quadratic form which corresponds to a semi-linear continuous Hermitian form 

N N 

and L{u) a linear continuous functional defined on H. Moreover, since 

N 

iiu)>Y, irl'\x))-'\u^\x)\'d^. 



^ f 2 2 

+ / {rf\x))~'^\uf'\x)\'^d'fi>c\\u\\H E c=const, 



we obtain, using Remark 1.1 to Theorem 1.1 from [I], that there exists one and only 
one element u = {u^i \ • • • , u^^\ uf\ • • • , u^^^) G H such that mp 



I{u) = inf I{u). 

Therefore, for any fixed v E H and r G M^, the function s(r) := 
one minimum point r = 0, so that 

d 

—I{U + TV) 1=0 = 0. 



I{u + Tv) has only 
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This yields 



= l^I{u + TV) l^^o = lim ^{liu + TV) - liu)) 



= lim^ / q ^(x) (z{x;u + Tv)z{x;u + Tv)—z{x;u)z{x;u)] dT 



-(1)/ \ 



^ N „ r 



-uf\x)uf\x) 



Calculate the first limit in the right-hand side of the last relationship. Taking into 
notice the notation for z{x]v) and the equality z{x]u + tv) = z{x]u) + tz{x;v), we 
have 

lim— / q~^{x) [z{x;u + Tv)z{x;u + Tv) — z{x;u)z{x;u)] dT 
r— »o 2t Jp V / 



1 



dr 



= lim— / q^\x) 2tRc{z{-,u)z{-,v)) + 0{t^) 
r^o 2r Jp L 

= Re J qi'^{x)z{x;u)z{x;v)dT. 

Performing similar calculations for the remaining limits we obtain 

r N 

= Re / q^^{x)z{x;u)z{x;v)dT + y2^^ / irl^\x))~'^u\^\x)vl'^\x)dji 

+ ^Re {rf\x))-''uf\x)vf\x)d-fi. (2.91) 

On the other side, for any fixed v ^ H and r G R , the function si(ri) := I{u + iriv) 
has the unique minimum point at ri = 0, so that 

d 



dri 



I{u + iTiv) I = 0, 



which yields 



n ^ r 1 1 — 

= Im / qi'^{x)z{x]u)z{x]v)dV + \^\m I {rf\x))^'^uf\x)vf\x)d'^. 
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^ f 2 2 ~ 

.■1 J^, 



From the latter, in line wil;Ii\ l2.911 ) and ( 12.921 ) . it follows 



0= / qY'^{x)z{x]u)z{x]v)dT + / {rf\x)) 

i=l "^T^ 



u\ {x)v\ {x)d'ji 



2.92) 



(2.93) 



Introduce function p{x) as the unique solution to the problem 

(A + k'^)p{x) = in Q', 
— = ^1 z[-]u) on r, 



dp{x) 



dp 
du 



0, i=l,N, 



dr 



ikp{x) = o(l/r), r = r oo. 



(2.94) 
(2.95) 
(2.96) 

(2.97) 
(2.98) 



Transform the first term in the right-hand side of (12.911) using equalities (12.941) — ( j"2. 981) 
and applying the second Green formula in domains Hj, i = 1, A^, and ^}R to functions 
p and z{-; v). We have 

= / —{A + k^)p{x)z{x;v)dx+\^ I —{A + k'^)p{x)z{x;v)dx 



i = l '^^i 



dp 



= / —{A-'rk'^)z{x;v)p{x)dx— I z{-;v)^-^ dT — T,ji{z{-;v),p) 

N „ N „ 

+ / —{A + k'^)z{x]v)p{x)dx — '^^ / 



dz{-]v] 



i=l 



=1 ^7. 



p{x) - z^{-]v) 



dp 



N 



dz{-v] 
dvA* 



p 



■dp{x) 
;v)— I dji 



dvA 



■dp 



/ —{A + k'^)z{x\v)p{x)dx— / z{-\v)^ dV — Tiji{z{-\v)^p) 
'n' JT dv 



N 

i=l -^^^ 



dp 
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N 



+E 



dz{-;v] 



dz{-]v] 



=1 ^7. 



dvA* } \ duA* 
= - j v)q{'^{x)z{x] u) dV - ^r{z{-; v),p) 



Calculating the limit as ^ oo and taking into account that Hr{z{-; v)^p) = o(l), we 
obtain 



,(1) 



(2,1), 



,(2)/ 



E / ^(^) / [K['^'\^.x)vt\o+Kr\^,x)vr\o]d^,^d^,^ 

tm?f [Af'^'(e,x);F¥ + A-P'K,-)?^]<i7,<i7. 

„• 1 i/Ti ^Ti 



Next, by virtue of i^M), 



N „ N „ 

^ / (r<"(,T))-^u|"(x),,<"(x)d7, + ^ / (rf (,T))-^ttp'(x)«f'(^)rf7 

! = 1 •'T 1 = 1 



TV 



z-l ^7. 



(2,2), 



Rewrite the last equality in the form 



+Af'^>(.x.f)^)<i7, 



(2.99) 
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+ (x) - (rf {x))W (A-f"{x,«pK) + 



Setting in 



i = 1, A^, we find 

N 



f2 f {rf\x))-' v!f>(x) - {rf>{x)f [ (/if '"(x.^pK) + 

i=l •^7i 



diy J 



N 



+ 



and consequently, 

y!t\x) = {Tf\x)f 

uf\x) = {Tf\x)f 



(1,2), ,.dp{i) 



dv 



C^7z.o 



^7^^, i = l,N. 



Substituting these quantities to ( 12.591 ) and ( 12.621 ), setting z{x) = z{x]u), and taking 
into account ( 12.941 ) — ( [2.981 ), we arrive at problem ( 12.731 ) — ( [2.821 ): the unique solvability 
of this problem follows from the fact that functional ( 12.641 ) has the unique minimum 
point u. 

Now let us establish the validity of ( 12.841 ). Substituting expressions ( 12.721 ) to ( 12.641 ), 
we obtain 

Cr^ = I{u^i\ . . . , u''^\ uf\ . . . , wj^^) = 
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= J qi\x)\z{x)f dr 
= J qiHx)\z{x)\'^ dr 

N 

N 



(1,2), 



d% 



(2,2). .MO 



Next, using relationships ( 12.721 ) — ( [2 ■821 ). we have 



d% 



di^^. (2.100) 







N r 

{/\ + k^)p{x)z{x) dx + '^^ I —{/\ + k^)p{x)z{x)dx 

i=i '^^i 



— r uuix) 

■{A + k'^)z{x)p{x) dx — / z{x)^- — 



dr - ^b{zjP) 



gA((9/-'-S)^' 

TV 



N „ TV „ 

+ E / -{A + k'^)z{x)p{x)dx + ^ 
1=1 -^^^ i=i -^^^ 



dz\ _ dp\ 



I —{A + k'^)z{x)p{x)dx — J z-^ dV — Tiji{z,p) 



N 



-E 



1 ^7. ^7. 



rdz\ _ rdz\ 



pd^. 



I Xuo{x)lo{x)p{x) dx - q^'^lzl'^ dV -J:r{z,p) 
Jn' Jt 



N 

E 



dp{x] 



^1^7. A. 



c^7*, c?7z« c?7*. 
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N 



X Ki'^'\^. y) + ^-^Kf^'\i. y) ^7., ^7., ^7.. 



TV 



i=l "'Tj Jit Jli 



X 



X Kf^ii. y) + ^-^Kf^'\i. y) ^7., ^7., ^7.. 



= / lo{x)p{x) dx — / ^l^p (ir — Eij(z, j9) 

N 



7« 



c^7* X 



X 



N 



It 



dv 



X 



0^7*. 



= / Iq{x)p{x) dx — I qi'^\z\^ dV — Yjji{z^p) 



N 

E / / 



(1,2). .x^MO 



d% 



dv 



d% 



d% 



(2.101) 



Equality ( 12.841 ) follows now from two relationships ( 12.1001 ) and ( 12.1011 ) and the fact that 
^r{z^p) = o{l) when R 00. □ 
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An alternative representation for the minimax estimate in terms of the solution to 
a system of integro-differential equations is given in the next theorem. This solution is 
independent of the specific form of functional (12.121) . 

Theorem 2.5. The minimax estimate of (i2.12\) has the form 

[{^) = l{(p)= j lo{x)ip{x)dx, (2.102) 

where function <-p{x) is determined from the solution to the problem ^2.1031) — ^2.11S\) : 

^GL2(S,Fi,(n',A))0, (2.103) 

(2.104) 

= on r, (2.105) 



{/\ + k^)p{x,uj) = {) in Vl\ 



du 



,T.u;))],. = - / Af''>K.,T) [(r<"K))Vf'K,'^)-«f'K,^) 

J -7,- 



d% 



i2)^c\\2.,(2) 



d% 



dp{x, to)) 




du 


-j 



-(1,1), 



-(2,l)/t r/^(2)/tNx2,,(2), 



+ / Kr>{i.x) {rr\OYyr{iM-vf\iM 



dp{x^ u) 



dr 



on 7^, i = l,N, 
+ ikp{x, (u;) = o(l/r), r=|a;|, r^oo, 



(A + /c2)(^(x,a;) = m Q', 

= Qi p{-,uj) + hQ on 1 , 



dv 

•,cj)]^^ = 0, 
d(p{x^ uS) 



dv 



0, % = \,N, 



dr 



ik(p{x,Lj) = o(l/r), r = |x|, r oo, 



13 



2.106) 

2.107) 
2.108) 
2.109) 
2.110) 

2.111) 
2.112) 



Here L'^{^, H^^ci^' , A)) denotes a class of functions the belong to 1^(1], A)) for any R>0. 
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where 

«f'K,-) = {'-f'{?))V 

J 7,- 



(1,2).. ^d(p{ri,uj] 



du 



(2,2).. .d(f{r],uj) 



dv 



dl^,, (2.113) 



^7^„ (2.114) 



awe? the right-hand sides in ({2.1061) are considered for every realization of random func- 
tions yl^\^) = yl^\^.,uj) and yi'^\^) = yf\^.,u) which belong with probability 1 to the 



space L {'ji), i = l,N. Problem (i 2. 1031) — (i2.112\) is uniquely solvable. 

Proof. The proof of this theorem is similar to the proof of Theorems 1.2 and 2.4. □ 



Remark 6. Function (fi{x, u) which is determined from the solution to problem ({2.1031) — ({2.11 
can be taken as a good estimate of the unknown solution (p{x) to the initial Neumann 
problem ({2.1\) — f27B) (see Remark 1 on p. 

Set 



+Kf^'\^.v)Kf^%J){rl'\Of] d7.„ (2.115) 



+A'p)(^,r/)A'p)(^,x)(rf)(0)'l ^7., (2-116) 



+Kl'^'\^,v)W\7)ir?\0f] ^7.., (2.117) 



Kf''\x,r]) = 



+Kl'^'\^.v)Kf^\^){rf\0)'] ^7z„ (2.118) 



Then Theorem 2.4 can be formulated as follows. 
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Theorem 2.6. The minimax estimate of the value of functional l{(p) has the form 

^ N „ 

^(^) = E / (4\^)yl'\^) + uf\x)yf\x)) dji + c, (2.119) 

where c = Jp ^5^0 o^r, 



vf\x) = [rf\x)f[ 



/MO 

dp 

dm 



(2.120) 



(i7ij, i = l,A^, 



one? functions z and p are determined from the uniquely solvable problem 



- (A + P)z(x) = Xuoi^Wx) in Q', 



n r 
— = U on 1 , 



(2.121) 

(2.122) 
(2.123) 



'7i 



(1,2), 



dv 



dj,^, onj„ i = l,N, (2.124) 



dz{x) 


-n 


du 





(2.125) 



— + ikz = o(l/r), r = Ixl, r — ^ 00, 
pGFii,(f^',A), 
(A + k^)v{x) = in n', 

— = Q zq onL, 
ov 



\pix)]^, = 
dp 



dp{x) 
dv 



= 0, on 7i, i = 1, A^, 



J 7i 



5r 



— i/cp = o(l/r), r = r — ^ 00. 



(2.126) 
(2.127) 
(2.128) 
(2.129) 

(2.130) 
(2.131) 
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2.4 Integral equation systems whose solutions are used to express minimax 
estimates 



In the previous section, we have obtained the integro-diffcrcntial equations whose so- 
lutions are used to express minimax estimates. In this section, we use the developed 
potential theory in Sobolev spaces and reduce these integro-differential equations to 
integral equations over an unclosed surface which is a union of the boundary of domain 



O and surfaces 7^, z = 1, A^^, on which observations are made. This reduction allows one 
to decrease the dimensionality of the problem of finding minimax estimates. 

We define first, in addition to the single- and double-layer potentials introduced in 
the previous sections, the corresponding boundary integral operators Sk, K^j -K"^, and 

iSk<f)ix) -.= 2 J ^k{x,y)ip{y)dTy, xeV, 



in the three-dimensional case, their kernels arc determined by the formulas 



{x-y.-Uy) 

[l - ik\x - y\) 

\x — y\ 



=M^.y) \^ itk\x-y\-l), 



(2.132) 
(2.133) 
(2.134) 
(2.135) 



\x - y[' 



(2.136) 
(2.137) 



d'^<^k{x,y) 
duxd^y 



= ^k{x,y) 



{l-ik\x-y\) 



\x - y\ 
{x-y,-Ux){x-y,-Uy) 



\x-y\ 



(—3 -I- 2>ik\x — y\-\- k'^\x — yp) 



in the two-dimensional case, 



d 



4\x — y\ 



4:\x — y\ 
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(2.138) 



(2.139) 
(2.140) 



d^M^.y) ikH^^\k\x - y\) I ^ 2 

r. r. = J. r. [Vx,Vy)\x-y\ 

ouxOVy 4|a: — yy^ L 

-[x - y, Ux)[x - y, iyy){\x - y\ + 1) H ^ {x - y, Ux)[x - y, Uy), 

J 4|x — yf 

(2.141) 

where H\{z) denotes the order-one Hankel function of the first kind and (•, •) the inner 
product in R", n = 2, 3. 

Note that, for example, in the three-dimensional case, the kernels of integral opera- 
tors ( 12.1321) — ( 12.1341) have a weak singularity and the integral in the right-hand side of 
( 12.1351) is understood as a Cauchy singular integral. 

Let us formulate the properties of the operators introduced above that are essential 
for the reduction of problem ( 12.1211) — ( 12.1311) to a system of surface integral equations. 

If r is a C^-surface, then the following operators are continuous at \s\ < 1/2: 

Kk : H^I^^'{T) ^ //3/2+^(r), 

Tk : H^/^+'{r) H-^/^+'{r) (2.142) 

(similar statements are valid for the operators K_^, ^^k^ T~k)- 

Also, operators Kk and K'_j^ acting, respectively, from H^I'^{T) to H^^'^iT) and from 
i7~^/^(r) to H^^/'^iT) are compact according to ( 12.1421) and the following equalities 
hold: 

{Skip. ^)L^iT) = {<P. ^-fc^)L^(r) G H-'/\r), ^ e H-^'\T), (2.143) 

where 

{f.g)L^{T) = j^fgdV =: g{f) 

for every / G H^I'^{T) and g G H-^I'^{T). 

Denote by (Vp?/^)^ (x) and (Wp?/')^ {x) the restriction on the domain \ Q of the 
single- and double-layer potentials ( 12.191) and ( 12.201) with a density ifj G i7^(r), s G R. 
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The traces on F = dfl of these functions and their derivatives satisfy the relations 
pp. 224, 225: 

{V^^^^)^ = ^Sk^^^ in W+\r), (2.144) 

i^ri^f = l{i^ + Kki^) in H^T), (2.145) 

|;(Vr'^)^ = -^(^-^i^) in H-^iT). (2.146) 

§^ i^^r^f = l^ki^ in (2.147) 
Similar relations are valid if we replace k by —k. 

Denote by D{^}) a countable set of positive values of wave number k with a limiting 
point at infinity such that the homogeneous internal Nuemann problem 

Ai; + = in (2.148) 

V = Oonr (2.149) 



has nontrivial solutions. Then ^ 

: Ai; + /c'l; = in J^, v = onT} , (2.150) 



r 



where N{I — K'^ denotes the null-space of / — K'^^. 
It is known that the solution V\{x) of the problem 



Ai;i(x) + A;V(a:) = in R^\n, (2.151) 
dvi 



du 

dvi 



fi on r, (2.152) 



dr 

can be represented as 



ikvi = o(l/r), r = |x|, r ^ oo if Imk > 0, (2.153) 



v^{x) = W^ipi{x) - V^fi{x), xeM.^\n, (2.154) 

where the function ipi := -uilr which is the trace of this solution on F can be determined 
directly or as a solution to the integral equation 

ipi{x) - Kk(pi{x) = -Skfi{x), X G r, (2.155) 
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when k ^ D{Q); generally (that is, for any A; 7^ 0, Im /c > 0) it can be found from the 
equation 

ipi{x) - Kk^i{x) - iT]Tk(pi{x) = -Skfi{x) - ir]{fi{x) + K'Ji{x)), x G T, (2.156) 

in which a number 77 G M, r/ 7^ 0, is chosen so that 

7/ReA;>0. (2.157) 

The existence of solutions to these integral equations at any /i G i7^^/^(r) follows 
from the unique solvability of BVP (12.1511) — (12.1531) . The solution to (12.1551) is unique 
because 

dimiV(/ - Kk) = dimiV(/ - K'^^) = dimiV(/ - K'_-^) = (2.158) 

due to (12.1501) and the Fredholm alternative; the uniqueness for (12.1561) is a consequence 
of the fact that the operator I-Kk- ir]Tk : H^/^{T) H-^I'^{T) defined by the left- 
hand side of (12.1561) is an isomorphism [HT| . 

A reasoning similar to that in [SI] — [53] enables us to prove that the solution V2{x) 
to the problem 

/\v2{x) + k^V2{x) ={)m^^\^, (2.159) 
dv2 



dv2 



= f2 on r, (2.160) 



dr 

can be represented as 



+ ikv2 = o{l / r) , r = |x|, r ^ 00, (2.161) 



V2{x) = V2(a^) - Vf ^2(2:), xeM.^\n, (2.162) 

where the function (p2 '■= V2\r which is a trace of this solution on T can be determined 
directly or as a solution to the integral equation 

ip2{x) - K_iif2{x) = -S_if2{x), x G r, (2.163) 

for k ^ DiQ)] generally, this function can be determined from the integral equation 
(for any A; 7^ 0, Im A; > 0) 

ip2{x) - K^np2{x) + ir]T_np2{x) 

= -S_-j2{x) + ir]{f2{x) + K'_-j2{x)), xer, (2.164) 
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where ij satisfies conditon ( 12.1571) . 

Indeed, the solution V2 G \ Q) of the Helmholtz equation ( 12.1591 ) that sat- 

isfies radiation conditions ( 12.1611) and the property ^ G admits the integral 

representation 

V2{x) = I ^%|^^2(2/) dVy - I y)^^dTy, X e (2.165) 

replacing in this formula by /2, we obtain equality ( 12.1621) for the solution V2{x) 
to problem ( |2J^ - ( I2J6T1 ) . 

Calculating the traces on T for both sides of ( 12.1621) and using relationships ( 12.1441) 
and ( 12.1451) with k replaced by —A;, obtain a boundary integral equation for (p2 

^2 - K_i^2 = -S^kh on T; (2.166) 
taking into account ( 12.1461) — ( 12.1471) . we can obtain another integral equation for ip2 



r_fc(^2 = /2 + i^V2 = on r. (2.167) 

Multiplying both sides of ( 12.1671) by a number rj and adding to ( 12.1671) . we find that 
ip2 satisfies integral equation ( 12.1641) . 

The uniqueness of solution of this integral equation follows from the fact that oper- 
ator K_j, is adjoint to K'j^ and operator T_i is adjoint to T^; therefore, / — K_j, + irjT_i 
is adjoint to / — K'j^ — iijTk. Also, 



Im (/ — K_j^ + irjT_i) = Ker {I — K'j^ — irjTk) 



_L 



In [S2] it is proved (see Theorem 3.34) that under condition ( 12.1571) Ker {I—K'f^—iijTk) = 
0. Thus, if ( 12.1571) holds, then integral equation ( 12.1641) is uniquely solvable. 

In order to reduce problem ( 12.1211) — ( 12.1311) in unbounded domain R^\n to a system 
of boundary integral equations let us apply the results formulated above in Subsection 
2.5. 

Introduce functions gf'\uf'\uf'') G HQ(^('ji) and gf'\u^^\uf^) G L'^i'Ji) defined on 

(1)/.(1) 42)x (1)/ .(1) .(2)x 



d7^,. (2.168] 
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Ki'''\^r)4\0 + Kr'i^:-)ur\0 ^7.,, (2.169) 



(2,1).. x.(2). 



Introduce also a function Zin = Zin{-',u) G Hl^^{W^ \ {[J^iji),A) which solves the 
BVP 

- (A + P)z,n = X.,{x)k in W \ (U£i7,), (2.170) 
dzi, 



du 



= pf\uf\uf^)on^,, i=l,N, (2.171) 



dzif^ 
dr 



+ ikzin = o(l/r), r = |x|, r ^ oo; 



(2.172) 



in line with Theorem 12 .21 . this function, in the domain \ (U^^Ti); determined 
according to 



+ /" $„fc(-,2/)/o(?/)% (2.173) 



Then the solution to problem (12.1211) — (12.1261) can be represented in Q! as 



Z Zg ~\~ Ziji^ 



(2.174) 



where the function z^ = Zs[-\u). 



u = (u 



-(1) -(2) 



(2), 



6 (i'(7l) 



X . . . X 



is determined from the solution to the following problem 



ZgeHU^'.A), 

(A + P)z, = in Q', 
dza dz 

on 1 , 



dzs 
dr 



dp du 
+ ikzs = o(l/r), r = r cxd, 



(2.175) 

(2.176) 
(2.177) 

(2.178) 
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where 



dv \ dv dv ] dv 



+ / ^^^^k{y)dy. (2.179) 



'Wo 



According to ( 12.1621 ) and ( 12.1641) . where V2 
trace of Zs on F denoted by 



Ir 



:= Zs\t = z\y - z, 



see ( 12.1741) ) satisfies the integral equation 

d Ziji 



Introduce also the notations 



du 



■ / dzifi 



and ip2 = ^s|r, the 



(2.180) 



Z = l,iV. 



(2.181) 



(2.182) 



Then, by virtue of (ESI) and ( ICTTD |^ = g^-^z = ^(^/^ + z,n) on T. 



Taking into account relationships ( 12.1541) , ( 12.1561) in which we set fi = Qi z = 



9i ^('^ + Zin), (fi = X, and vi = we obtain an integral representation for the solution 
p of BVP I^M)-I^M) in the domain R^\Q 

p = - Vkl\i^ + Ir), ^ = MV, 



where X = J^lr satisfies on F a boundary integral equation 

{I-Kk- irjTk)x = -{Sk + ir]{I + K))qf{i; + z 



(2.183) 



(2.184) 



From ( 12.1831) it follows that 



dp 
dv 



du 



dV^q^'^{ilj + z^n \y] 



du 



, i = l,N. 



and 



pL = L, - VrV(^ + Ir) L, , ^ = 1, A^, 



«n ir; I7, 



(2.185) 



(2.186) 
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Since (p'P = and (pf^ ~ ^ ? ^ = 1? ^? the latter equalities combined with ( 12.1811) 



and (EIHD mean that functions t/', x, and \ defined by (EHUD, (ElMl) 

and ( 12.721 ) solve the following integral equation system 

dZr 



dv 



{I - Kk- ir]Tk)x = -{Sk + ir]{I + Kj.))q^ + z,n lr) 
= Wr'x I,. - V^C'(V^ + lr) I,. , ^ = T;A^, 



(2) dW^x 
(/}■ = — 

dv 



dp 



(1) _ /.(l)^2 



w; = r 



kf"(-.«)^!"K) + Af'^'(..0^f'K)]d7,. 

J 



(2.187) 

(2.188) 
(2.189) 

(2.190) 
(2.191) 



(2) _ /^(2)n2 



It; = r; 



'7i 



i^P^(-,o^r(o+i^r^(-,o^r(o ^t.,, ^ = i,a^. (2.192) 



-(2,2)/ tN, (2) 



Resolve this system with respect to functions u'P and uf\ Replace in ( 12.1681) and 
(E369]) uf^ and uf^ by their expressions ( |2J^ and ([2392]) to obtain ^f^(wf\wf^) = 
/^S^^l'/^S^Wr) and ^f^(wf\wf^) = pf\^f\^f^) where 



-(1,1), 



,(1)/ 



(^'')(.,r/)^f(r/)J (i7.„ 
i^P(.,r/)^W(r/)+i^p)(.,r/)^f(r/)l ^7.,, 



(2.193) 



(2.194) 



and k'^^'^\-^ = 1, 2, are determined according to ( 12.1151) — ( 12.1181) . Next, replacing 

in (12173]) functions g^\4\ uf ) by p^\^l'\ and by (^{^, ^f^) 

we have 



2,- 



^fc (2), (1) ^ (2) 



fc fl). fl) , (2). 



(2.195) 



Substituting this expression into ( 12.1871) — ( 12.1901) we conclude that il)^ x-, '^T' 1 and Lpf^ 
satisfy the integral equation system 



Si: 



du 



+ 



dN_kk 
du 



(2.196) 
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{I-Kk- ir]Tk)x = -{Sk + iriil + K)) 



+ iV_fc/o|r , (2.197) 



(1) 



r 



, i = l,N, (2.198) 



(2) ^ ^ 
dv 



i = l,iV. (2.199) 



Summing up the above reasoning and taking into consideration Theorems 12.41 and 
2.61 we arrive at the following result. 



Theorem 2.7. The minimax estimate of l{(p) has the form 



N „ 

^) = J2 (^^\^)y\'\^) + uf\x)yPix)) + c, (2.200) 



where 



c = 



N 



i=l 



,(1)/ 



-(1,2), 



go (ir, 
(2.201) 

(2.202) 



uf\x) = {rl"'{x)) 



(2)/^XX2 



The auxiliary function 



iV_^/olr 
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defined on T and functions x •= J^lr, ^f^^ pl-y,, and (f'p := i = l^N 
determined from the solution to the integral equation system (12.1961) — (12.1991) . in which 
J] is an arbitrary real number such that r/Re /c > 0. This system is uniquely solvable for 
all values of wave numbers k, \mk > 0. 
The estimation error u = l^vfl^, wh ere 

- W-^i' V!", ') ^ + N_-,lo \r) (2.203) 

in the domain \ Q. 



are 



r 



Proof. It is necessary to prove only the unique solvability of system (12.1961) — (12.1991) 
which follows from the unique solvability of the system of integro-differential equations 
(|Z72])-(|Z83]). 

Indeed, let the integral equation system (12.1961) — (12.1991) has another solution i/j, 
X, (pf \ i = 1,A^. Introduce functions u\^\x), uf\x), p, and Zin by formulas 
(EHH), (E1921), (E203]), and I^JM) in which x, and ipf\ i = MV are re- 
placed by tjj, X? and (ff'\ i = 1,N, and functions Zg and z by Zg = Wp^'i/i — 
Vp^5m|r and 5 = + Zin- Then from (12.1981) — (12.1991) it follows that (f^^^ = p\j^ and 
(ff^ =1^ , i = l,N. Theorem [22] and the equalities gl^\ul^\uf^) = p\^\(fl^\ (ff'^), 

Qf\uf\up) = pf\(pf\(pf^) imply that z and p will also satisfy integro-differential 
equation system (I2.72l) - p.83l) which is uniquely solvable. The latter statement and 
the fact that operators I - Kk - ir]Tk, I - K_i + iT]T_i : H^I'^{T) H-^/'^{T) are 
isomorphic mappings yield (pf^ = ^f'P-, (p'P = <^f\ i = l^N, and ijj = ijj ^ x = x- D 

Using Theorem 12.51 the notations 

N 

4^ := f A-P'«.-)(rf'(?))Vf'K,^)d7.+ / A'F'K, ^(rf '(f))Vf'«,^) ^7.. 
7 - I n I ~ N ~(1) N ~(2) 5^(0 
that led to the proof of Theorem 12.71 we can prove the following 
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, z = 1, A/^, and the reasoning 



Theorem 2.8. The minimax estimate of I (if) has the form 



l{ip) = l{if), 



(2.204) 



where 



N 



i=l 



+ L>_fc Ir ) + /io 



(2.205) 



in the domain M^\Q and functions tjj, o-'^d '~pf \ = 1? ^5 o-f^ determined from 

the solution to the uniquely solvable integral equation system 



Sill 



+ 



dD 



-k 



dv 



(2.206) 
(2.207) 



N 



i=l 



fc n).~(i) ~(2). 



TV 



fc fl),~(l) ~(2), 



i = l,N, (2.208) 



~(2) WfcX 



7« 



-fc f2).~(l) ~(2) 



fc a).~(l) ~(2)^ 



i = l,iV, (2.209) 



where rj is an arbitrary real number such that rjKek > and functions p^p{(pf\(pp) 
and pf\(p^^\(pf^) are determined by (12.1931) and (12.1941) in which (p^^ and ipf ^ should 
be replaced, respectively, by and (^^^ 
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Note that integral equation systems ( I2.196l )-f l^l^ and ( I2.2()6I )-( |^:^IM ) are singu- 
lar. 

U k ^ D{^}), then repeating the reasoning used in the proof of Theorems [2/71 and 
Em where equalities of the form ( 12.1561 ) and ( 12.1641 ) are replaced, respectively, by those 
of the form ( 12.1551 ) and ( 12.1631 ) we see that the minimax estimate of /((/?) may be found 

(1) ,.(2) 



from (E^nni)- (122021) or (122D1, (122nSD, where the functions i/j, x, ^i , ^1 , ^ = 1, A^, 
and -0, X, \ i = 1,N, are determined from the solutions of the weakly singular 

integral equation systems ( I2J961 ) - ( |2:T991 ) and (ESQSD - (E^QSl) with r/ = 0. 



Remark 7. The assumption that surfaces 7^, i = 1, N, are pairwise non-overlapping 
is not essential. Slightly changing the proof, one can extend all results of this chapter 
to the case when surfaces "fi intersect on a finite system of contours. 

Remark 8. The method proposed in this chapter enables one to solve the problem of the 
minimax estimation of the value of a functional defined on ^{x,t) := Re [e^*'^V(x)] 
of the form 

L($) := / / lo{x,t)^{x,t) dxdt 

from the observations 

yf\x.t) = I K^'''\x,m^,t)d^,^+ [ Ki'^'\x,0^^^dj,^ + rjl'\x,t), 

yf\x.t) = I Ki'^'\x,ym,t)d^,^+ f Ki'^'\x,y)^^^dj,^ + rjf\x,t), 
Jii Jii 

i = 1, in a time interval from t = to to t = T. Here we assume that for f G Gg, 
a; > 0, /q £ L'^i^o x (to^T)) is a given function, and T]f^{x^t) and r]f\x.,t) are 
observations errors which are realizations of random fields defined on 7^ x (ig, T) that 
are continuous in the mean-square sense and have zero expectation and unknown second 
moments 'Ei\t]'P {x^t)\^ and'Ei\rjf'\x^t)\^ satisfying the inequality 



"^^0 Jit 



+ E / / E|^f^(:^,^)l' (r?\x.t)) d-f,dt < I, 



where rf\x.,t).,rf\x.,t) are given functions continuous on 7^ x (to^T), i = 1,N, that 
do not vanish on these sets. 
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PART 3 



Minimax estimation of the solutions to the boundary value prob- 
lems from point observations 

In the previous chapters we looked for estimates of unknown solutions (and the right- 
hand sides of equations entering the statements of the corresponding problems) from 
the observations of these solutions distributed on a system of subdomains or surfaces. 
In this chapter, we consider similar problems in the case of point observations and 
propose constructive minimax estimation methods. 



Let x'^, k = 1,N and Xk, k = l,m he given systems of points belonging to domain 
\ Q. The problem is as follows: to estimate the expression 

m 

K<p) = ^aiif{x^), (3.1 

i=l 

from the observations of the form 



Vk = ^{x'k) +r]k, k= 1, N, (3.2) 



that correspond to the system state (p described by problem (12.11 ) — ( l2~il ) in the class of 
estimates 



N 



K^) = '^u^yk + c, (3.3) 

fc=i 

linear with respect to observations (13.21 ) under the following assumptions: h G Gq and 
r] := (?7i, . . . ,?7iv) G Ci, where the set Go is given by formula (12.81 ). iji are errors of 
observations (13.21 ) that are realizations of random quantities iji = r]i{uj), G\ is the set 
of random vectors r/ = (r/i, . . . , r/^v) with zero expectations and finite second moments 
satisfying the condition 

N 

i=l 

and G C, ? = 1,A^,, G C, ? = 1, m, G M, z = 1,A^, and 7^ are given 
numbers. 

Set u := (wi, . . . , un) G M^. 

Definition. The estimate 

N 



i=l 
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in which numbers iii and c are determined from the condition 



sup E\l{(f) - ^ inf , (3.5) 



where 



N 



K'f) = ^u^y^ + c, (3.6) 



m = ^{xi) + r^z, 1 = 1, N, (3.7) 

and (p{x) is the solution to the Neumann BVP at h = will be called the minimax 
estimate of expression ^S.l\) . 
The quantity 

g:={ sup E|/((^)-/(^)p}i/2 (3.8) 

will be called the error of the minimax estimation ofl{(p). 

Based on the proof similar to that of LemmaEH] (in fact, much simpler) we can show 
that the following statement is valid in the case of point observations. 

Lemma 3.1. Finding the minimax estimate of functional l{(p) is equivalent to the 
problem of optimal control of the system described by BVP 

m N 

Az{x] u) + k^z{x; u) = ai5{x — Xi) — UkS{x — x'l^) in \ 

i=l k=l 

dz(-;u) 

= U on I , 



dv 

dz{-; u) 



-\- ikz{-;u) = o{l/r), r = \x\^ r ^ oo. 



dr 

with the cost function 

N 



/(.)= /,rW{.;«)rfr + V,--KI^-min. (3.9) 

1=1 

Starting from this lemma and proceeding with the reasoning that led from Lemma 
27n to Theorems 12.41 and 12.51 we arrive at the following result 
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Theorem 3.1. The minimax estimate of I (if) has the form 

K^) = ^Uiyi^c= l{(f), (3.10) 

where 

Uk = r\p{xj^^ /c = l,A^, c = J 'zhodr, (3-11) 

the functions z^p, G V'{M^ \ Cl) and (p = (p{-.,u) G T>'{E^ \ are determined, respec- 
tively, from the solution to the following problems: 

m N 

— (A + k^)z{x) = ai6{x — Xi) — r\p{x'^)5{x — x'f^) in Q, (3.12) 

i=l k=l 

Oz 

— = Oonr, (3.13) 
oiy 

Oz — 

— — \- ikz = o(l/r)^ r = \x\, r ^ oo, (3-14) 
or 

Ap{x) + k'^p{x) = in (3.15) 

^ = q^^z, on r, (3.16) 

dp 

— ikp = o(l/r), r = r oo. (3-17) 



and 



dr 

N 



- (A + k^)p{x) = J2 rlivK) - A4M^ - 4) ^ri (3.18) 

^ = Oonr, (3.19) 
ou 

dp 

— — \- ikp = o(l / r) , r = \x\, r ^ oo, (3.20) 
or 

A(p{x) + k'^(f{x) = in (3.21) 

^ = qfp^ho onT, (3.22) 

— ik(p = o(l/r), r = Ixl, r ^ oo. (3.23) 

Or 

Problem (iS.l^) — (iS. 23\) is uniquely solvable. The following estimate is valid for the 
error a of the minimax estimation of l{(p) 

/ m \ 1/2 

a=[l{p)]'/'=i^^-a,p{x,)j . (3.24) 
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In conclusion, we formulate the statements similar to Theorems 12.71 and 12.81 that 
enable one to reduce, in line with the algorithm applied in the proof of Theorem I3.ll 
the determination of minimax estimates to a problem of less dimensionality. 



Theorem 3.2. The minimax estimate of l{(p) has the form 



N 



(3.25) 



where 



k=i 



N 



-k\ 



1=1 



(3.26) 

g,dT, (3.27) 



and functions 

i/j := z\j. 



N 



,7 = 1 



00 q 



^qh{A)^--k{- 



1=1 



X '■= p\r, and numbers p{x'i) , I = 1, N, are determined from the solution of the following 
equation system: 



m 

E 



du 



N 



r 1=1 



{I - Kk- iT]Tk)x 



m 



N 



-Sk - iriil + Kl))q^^ ( ^ + ^-fc(- " " Y.rMx[) 



X 



1=1 



(3.28) 
(3.29) 

(3.30) 



m 



N 



in which rj is an arbitrary real number such that rjRek > 0. This system is uniquely 
solvable for all values of wave numbers k, Im/c > 0. 
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Equation ( 13.301) may be rewritten in the form 



N 



[1 - a^i\p{x'i} + ^ aup{x'i) = WkXi^d - H^i ^(a^z) + A, ^ = 1, (3.31) 
where 



1 2 / e 



zA:|4-y| 



27r Jr - ?/l 



-ifc|y-a::;| 



\y - x[ 



1 



ifc|x',-y| - 



-ik\y—Xj 



27T ^ \lr\xs-y\'^' '''' \y-Xj^ 



Theorem 3.3. The minimax estimate of I {if) has the form 

l{cp) = l{(p) = ^a,(^(x,), 

where 



dVy, s = l,N. 



(3.32) 



N 



and functions 



qlHi^ + - ^(^!)) - x'i)\t + ^0 ^n R3 \ Q (3.33) 



r ' 



1=1 



X '■= ^\t, and numbers ^{x'l), I = 1,N, are determined from the solution of the follow- 
ing equation system: 



{I - K'r + ir^T_-,)i 



N 



5$. 



00 1 



1=1 



dv 



{I - Kk- ir]Tk)x 



, (3.34) 
(3.35) 



= {-Sk-zrj{I + Kl)) 



N 



-k\ 



\ 1=1 



X 



(3.36) 



N 



in which ij is an arbitrary real number such that ijRek > 0. This system is uniquely 
solvable for all values of wave numbers k, Im/c > 0. 
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Remark 9. Similar results can be obtained if the estimated functional has the form 
(EH). 
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